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Abstract: In this paper we analyze perturbatively a gcjfi classical field theory with and 
without temperature. In order to do that, we make use of a path-integral approach de- 
veloped some time ago for classical theories. It turns out that the diagrams appearing at 
the classical level are many more than at the quantum level due to the presence of extra 
auxiliary fields in the classical formalism. We shall show that several of those diagrams 
cancel against each other due to a universal supersymmetry present in the classical path in- 
tegral mentioned above. The same supersymmetry allows the introduction of super-fields 
and super-diagrams which considerably simplify the calculations and make the classical 
perturbative calculations almost "identical" formally to the quantum ones. Using the 
super-diagrams technique we develop the classical perturbation theory up to third order. 
We conclude the paper with a perturbative check of the fluctuation-dissipation theorem. 
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1. Introduction 

Before starting this paper we would like to apologize with the authors of the famous paper 
[1] who first made use of the ancient-italian word " diagrammar 1 '' for having copied their 
idea. We used that same word because we feel "diagrammar" express in the most complete 
way the skills one has to develop to draw all possible Feynman diagrams of a classical 
theory which are many more than the quantum counterpart. We also think that the word 
express in a nice way the feeling we had when we realized that the very many classical 
diagrams could be encapsulated in few super- diagrams. 

Let us go back to the paper. In the last fifteen years the physics of heavy ion collisions 
has become one of the hottest topics in high energy physics [2]. The most interesting point 
is that this may lead to the formation of state of matter known as quark-gluon plasma 
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(QGP)(for a review see ref. [3]). This topic has attracted the attention of many sectors 
of theoretical physics, last but not least, even of strings using the duality AdS/QCD (for 
example see ref. [4] and references therein). 

We know that in the first instants after the (QGP) formation, the gluon occupation 
numbers is going to be very high [5], so one could describe the system by a classical 
field theory [6], analogously to what happens to QED which, when the photon occupation 
number is high, can be described to a good degree of accuracy by the Maxwell equations. A 
lot of work started in the nineties in this direction [7, 8, 9, 10]. Most of these authors studied 
classical field theories by just solving the Hamilton's equations of motion and developing 
from there perturbative calculations. Only the author of ref. [10] used a different route 
making use of a path-integral approach for classical mechanics developed in the eighties 
[11]. This is basically the functional counterpart of the operatorial approach to classical 
mechanics developed in the thirties by Koopman and von Neumann [12]. We know that 
the path-integral is the most natural tool to use in order to develop the perturbative 
calculations and the associated Feynman diagrams. It was so for the quantum path- integral 
(QPI) [13] and we will show that it will be the same for classical systems via the classical 
path-integral (CPI). In this paper, like in ref. [7] and Ref. [10], we will limit ourselves to 
a gcp 4 theory. 

The paper is organized as follows. In sec. 2 we will briefly review the CPI in order to 
make this paper self-contained. We shall show that, besides the phase-space variables, this 
formalism makes use of extra auxiliary fields. These have a well know physical meaning, 
being associated with the Jacobi fields [14] of the theory plus its symplectic conjugates 
and the response fields [16]. The presence of all these fields make the diagrammatic in the 
CPI rather complicated. In sec. 3 we develop the formalism in the case of a field theory 
without temperature and present the associated diagrams which, for a point particle, have 
already been worked out in Ref. [17]. In sec. 4 we develop the formalism in the presence 
of temperature and make contact with the calculations developed in Ref. [7] and Ref. [10]. 
We will notice that there are cancelations among various diagrams. This cancelation is 
due to some hidden symmetries present in the CPI. One of these is a supersymmetry and 
this leads to the use of super-fields (for a review see Ref. [18]) and their associated super- 
diagrams which we will introduce in sec. 5. In the same section, we will show how the 
super-diagrams will allow us to reproduce the results of sec. 4, but with fewer diagrams. 
These super-diagrams have the same vertices as in the quantum case. This is a consequence 
of the relation between the CPI and the QPI which was studied in Ref. [19]. In sec. 6 we 
conclude the paper with a perturbative analysis of the fluctuation-dissipation theorem [16] 
using the super-diagrams developed in sec. 5. We confine the detailed derivation of several 
results to few appendices. 

2. The Path-Integral for Classical Mechanics (CPI). 

In the thirties Koopman and von Neumann [12] developed an Hilbert space and operatorial 
approach to classical statistical mechanics. Their formalism is based on four postulates: 

I. A state is given by an element \ip) of an Hilbert space. 
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II. On this Hilbert space the operators q and p associated to the classical variables q 
and p commutes, i.e. [q,p] = 0. If we indicate with (p a = (q,p), a = 1,2, then the 
simultaneous diagonalization of q and p can be expressed as 



0>8> = ¥>>8>- 

III. The evolution of ij)(q,p) is given by the Liouville equation 

i ^ =Hip(q,p,t), (2.1) 

where % = —id p H(q,p)d q + id q H(q,p)d p is the Hamiltonian operator. 

IV. The Liouville probability density p(q,p) is given by 

p(q,p) = mq,P)\ 2 - (2-2) 

As a consequence of the four postulates above also p satisfies the Liouville equation as 
it should be. 

As this is an operatorial formalism there should be an associated path-integral and, 
in fact, it is easy to build one [11]. Let us start from the classical transition amplitude 
(tp a , t\(fQ, to) which clearly is 

(vA t ) = ~5[<p a - ^(t; cp , t )} (2.3) 

where <j) c \ is the solution at time t of the Hamiltonian equation of motion 

r) FT 

r-^e? < 2 - 4 > 

with initial configuration ipo at time to- Li Eq. (2.4) uo ab is the symplectic matrix [25] and 
H the Hamiltonian of the system. 

Slicing the interval of time in Eq. (2.3) in N small intervals and replacing the functional 
Dirac delta on the r.h.s of Eq. (2.3) with the expression 

5& a - c(> a ci(t; <Po, t )] = 5& a - cj ab d b H]det(5^d t - cj ad d d d b H), (2.5) 

where det[ ] is a properly regularized functional determinant, we get that Eq. (2.3) can be 
written as 

( V a ,t\^,t ) = [ £}"(p a £>\ a £} C a £>c a exp i [ drC 
J I Jt 

where @"<p indicates the integration over all points, except the end points which are fixed, 

while the symbol 3> indicates that the end points are integrated over. The £ is 

C = A a a + ic a c a - X a oj ah d b H - ic a u ad d d d b Hc b . (2.7) 

The auxiliar variables A a have been introduced to do the Fourier transform of the Dirac 
delta on the r.h.s of Eq. (2.5), while the Grassmanninan variables c a and c a are needed to 
exponentiate the determinant in Eq. (2.5). 



(2.6) 
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In Ref. [11] it was erroneously written that the determinant above is one. This is 
not so, but it depends on the regularization used to calculate it. If we discretize the time 
and do the same for the matrix associated to the determinant, then this turns out to 
be one, or a constant, if we use the the ltd calculus [20] or pre-point discretization [21], 
while it turns out to be non-constant and dependent on the ghost fields if we use the 
Stratonovich calculus [22] or mid-point discretization. It is only with this discretization 
that the Feynman rules which we derive from the continuum theory are going to be the 
same as those in the discretized case. This was shown in details in [23]. If we had used 
the Ito calculus in which the determinant is constant and could be disregarded, then the 
Feynman rules derived from the continuous version of the non-Grassmannian part could, 
in principle, not be the same as those derivable from the discretized version. This is the 
reason we use the Stratonovich prescription in which discretized and continuum Feynman 
rules coincide. In this prescription, the determinant depends on the fields and so there is 
a true interaction between the Grassmannian variables c a , c a and cp a . The vertex coming 
from this interaction has to be considered in order to get the right results. This is not 
the only reason to keep the Grassmannian variables c a and c a . Physically, they are the 
Jacobi fields and their correlations give the Lyapunov exponents [24] . These quantities are 
related to the matrix elements of the determinant above and to get these elements we need 
to couple c a , c a with external currents as explained in ref. [11]. 

Let us now go back to Eq. (2.6) and build the equivalent operatorial formalism. Given 
two variables 0\ and O2 from the path-integral (2.6) we can derive the commutators defined 
as follows [13, 14] 



Oi,0 2 



lim(Oi(t - e)0 2 (t) - 2 (t - e)Oi(t)). (2i 



The symbol () means the average under the path-integral and sandwiched among any state. 
We have put a hat " on the variables on the l.h.s. of Eq. (2.8) because they become operator 
via this definition. Applying (2.8) to the basic variables (p a , X a , c a , c a we get: 



<P a ,X b 



i5 a b , [c a ,B b }.=5 a b , (2.9) 



where by [ , ] + we indicate the anti-commutators. All the other commutators and anti- 
commutators are zero in Eq. (2.9). In particular, [i^ a ,(^ b ] = 0, which implies [q,p\ = 0, 
and this confirms that we are doing classical mechanics. From the first equation in (2.9) 
we get that a representation for A a is 

A* = (2-10) 

and this identifies A a with the response field (see for example Ref. [16]). 
Let us first write down the Hamiltonian associated to the Lagrangian C 

H = \ a u ab d b H + ic a uj ad d d d b Hc b . (2.11) 

Having obtained our operatorial formalism let us see if we can recover the KvN approach. 
To do that, let us consider the first piece of which we call Bosonic one (B) 

H B = X a u ab d b H. 
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Using Eq. (2.10) we get that % goes into the following operator 

= —idpHd q + id q Hdp = L 

and so we obtain the Liouville operator. This confirms that, at least for the Bosonic part, 
the path integral behind the KvN formalism is exactly the one presented in (2.6). As 
we said after the KvN postulates, both the elements of the KvN Hilbert space 4>{^>) and 
the probability densities p(ip) evolve via the same operator %. As a consequence also the 
path-integral for the evolution of both ip and p is the same. In this paper, as was done 
in the first papers [11], we will consider the averages taken with the probabilities densities 
p(<p) . The reason for this choice is due to the fact, as it has been proven in Ref. [15], 
the Hilbert space of KvN is actually made of just Dirac deltas on the phase-space points: 
5(ip — ipo) = (ip\ipo), and the superposition principle does not act on these states. This is 
natural in CM. So saying that we use a Hilbert space formalism in CM is a rather "formal 
"statement, and for this reason we prefer to use the probability densities. 

Going now back to T~L the curious reader may ask what is the mathematical meaning of 
T~L if we keep also the Grassmannian variables as in (2.11). Actually it has been proved in 
ref. [11] that it becomes a generalization of the Liouville operator know in the literature [25] 
as the Lie derivative of the Hamiltonian flow. This is easily proved once one realizes that 
the c a behaves under symplectic transformation as basis of differential forms of ip a , while 
the c a are the symplectic duals [25] to c a or, in other words, the basis of the totally anti- 
symmetric tensors. All the geometrical aspects of the path-integral for classical mechanics 
have been analyzed in details in [11, 26] . In the second reference it has been proved that 
the set of variables (p a , c a , c a , X a "coordinatize "a double-bundle over the phase-space [25]. 
The whole Cartan calculus [25, 26] can be reproduced using these variables and a set of 
seven charges present in this formalism [11]. Among those charges we will mention just 
two which are important for the rest of the paper 

Q H = ic a X a - c a 8 a H 

Q H = it a 0J ah \ h + c a u ab d b H. 

They make up a universal N = 2 supersymmetry for any Hamiltonian system. In fact, they 
close on the Hamiltonian 

Qh,Qh\=2iH. (2.14) 

We will see that, thanks to this supersymmetry, there will be many cancelations among 
Feynman diagrams once we embark on the perturbation theory. The geometrical meaning 
of these super-charges has also been investigated: they are basically the exterior derivative 
on the constant-energy surfaces of the phase-space. Or, in a more mathematically technical 
language, they are associated to the equivariant cohomology of the Hamiltonian field [27]. 

In this paper we will make use of another tool typical of supersymmetry which is the 
one of super-fields (for a review see ref. [18]). This is an object which allow us to put in 
the same multiplet all the variables ip a , X a , c a and c a . To do that we have first to extend 
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time t to two Grassmannian partners of time 9, 9. The super-field & a (t,9,9) is a function 
of t, 9, 9 denned as 

$ a (t, 9, 9) = ip a + 9c a + 9oj ab c b + i99uo ab \ b . (2.15) 

It is then possible to give a much simpler expression for the Lagrangian in C (2.7). In fact, 
it is easy [19] to prove that 



d9d9L[<$> a ] + (s.t), (2.16) 



where L is the standard Lagrangian associated to the Hamiltonian H of classical mechanics 
which we have used in the r.h.s of Eq. (2.7) and {s.t.) is a surface term. 

Using the super-field it is also possible to give a generalization of expression (2.6), that 
is the following (for its derivation see Ref. [19]), 

($ o ,t|$g,t ) = J @"<S> q m p expji J^idrd9d9L[^ . (2.17) 

The indeces p and q indicate the first and second of the indices of ip a . They are called 
this way because we know that ip 1 = q and ip 2 = p. 

An expression similar to (2.17) can also be written for the generating functional, which 
instead of being the complicated object: 

Z[J = 0] = J @<p a @\ a @c a @c a exp^i J t dT ^j ( 2 - 18 ) 

can be written as: 

Z[J = 0] = J 9)§ a exp ji dr J id9d9L [$] + (s.t.) j . (2.19) 

Note how similar is (2.19) to the quantum analogue: 

Zqm [ J = 0] = J &ip a exp 1 1 J drL [ip] j . (2.20) 

Basically ip a is replaced by <3? a and 1/h by i f d9d9 while the functional weight is given in 
both cases by the standard Lagrangian L of classical mechanics. 

The similarity between Eq. (2.19) and Eq. (2.20) has been studied in great details in 
Ref. [19]. In this paper this similarity will be used to drastically simplify the perturbation 
theory which, instead of making use of the complicated Lagrangian C of (2.7) appearing in 
Eq. (2.18), will make use of just L appearing in Eq. (2.19) borrowing many results from 
the quantum analog (2.20). 

Before closing this section we would like to bring to the attention of the reader the 
physics behind the variables c a . We do that because people may consider them useless. 
We have already made hints before to their mathematical meaning, being that of basis for 
differential forms (more details can be found in ref. [26]). Their physical meaning has been 
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explored in Ref. [24] and it goes as follows. Let us look at the equation of motion for c a 
which can be derived from the C of Eq. (2.7) 



This equation is the same as the one satisfied by the first variations Sip a : 

These first variations are also known as Jacobi fields. (See for example Ref. [14].) So we 
can identify c a ~ 8ip a . The correlations of the Jacobi fields are related to the Lyapunov 
exponents of chaotic systems: 



lim (6<p(0)6<p(t)) ~ e 

t— >-oo 



-At 



where (A) is the highest Lyapunov exponent of the system. This is equivalent to saying 
that the correlations of c a , c a gives the Lyapunov exponents [24] So we see that the c a , c a are 
crucial ingredients to get a full information on the dynamical system [24] we are studying. 
We talk about correlations between c a and c a because the correlations among only c a gives 
zero [33] [24] [34]. The c a are just their symplectic dual and, once they are multiplied by the 
symplectic matrix, they can also be identified with the Jacobi fields. So the correlations 
among c a ,c a are also equivalent to correlations among jacobi fields. Of course, the Lya- 
punov exponents could be obtained in many other ways [28] without using the c a ,c a , but 
we find this manner rather elegant. Actually, we know that classical mechanics could be 
built by just using the (p a , and so the variables A B , c a and c a are redundant variables and 
this redundancy is signaled by the symmetries present in C like the supersymmetry and 
other invariances. But, as usual in physics, the redundancies and the symmetries make 
things more elegant and allows the use of all the tools that group theory put at our dis- 
posal. In our case, if we had not used all the variables, we would not be able to build the 
super-field and use the expression (2.19) for the perturbation theory. 

Before starting the next section we should introduce for completeness a new represen- 
tation of the commutator algebra (2.9). More details can be found in Ref. [19]. This new 
representation will be very important in order to do perturbation theory using super-fields. 
We had seen that C in Eq. (2.16) and L[<J>] differ by a surface term (s.t) and the same 
will happen at the level of generating functionals (2.17) and (2.19). In order to get rid of 
that surface term the trick [19] is to change the representation in (2.10) that we derived 
from the commutator (2.9). We could, for example, represent q and X p as multiplicative 
operators and p and X q as derivative ones 



Analogously we can proceed for the Grassmannian variables: represent c q and c p as mul- 
tiplicative operators and c p and c q as derivative ones: 

<r = i± 

aCp (2.24) 
& d_ 

The generalized states [19] in the extended Hilbert space, including the forms, will then 
be (q,Xp,c q ,c p \. Note that the variables q, X p , c q , c p are exactly those which enter the 
super-field <fr q 

$v = q + 0c q + 9c p + ieeX p . (2.25) 

This is crucial in order to give a complete super-field representation of our theory without 
any surface terms. We invite the reader to master the details contained in ref.[19] before 
embarking in the rest of this paper. 

3. Perturbation Theory Without Temperature. 

What we will do in this section for a gcf) 4 scalar field theory, has been partly been done for 
a point particle in Ref. [17]. We use the word "partly" because those authors neglected 
the Grassmannian variables c a and c a . So their analysis apply only to the zero- form sector 
of the theory. We called this sector this way because, as indicated before, the c a are the 
basis of the differential forms [25], [26]. 

Let us start with a scalar field theory whose Hamiltonian is 



H = I d 4 x 



ir (V(p) 2 m 2 4> 2 <f 

- + v ' H — + q—r 

2 2 2 y 4 



(3.1) 



where tt(x) is the momentum conjugate to the field 4>(x). We will indicate with Jtff the 
Hamiltonian density: 

tt(x) [V0(x)] 2 mV(x) 6\x) 

^w = — + — ^ — + — 2 — g ^r ^ ' 

Let us next build the generating functional associated to the classical path-integral of the 
previous section, but we shall choose only those trajectories that start from a fixed point 
in phase-space ipi = (</>j,7Tj) and we will not average over this initial configuration. The 
expression is: 

2<fii[J<j>l J\*i Jon] = 

= J 9'(p®\9c9ceyLv^iS + i J J dt J d 3 x + Jx^K - iJ c4 ,c' l> - ic n J 5 ^ j ^ ^ 
where 

S = J J dt j d 3 x (\ a ip a + ic a c a - \ a oJ ab d b H - ic a u} ad d d d b Hc b ^j (3.4) 
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with the H given by Eq. (3.1) and the symbol & indicates the "sum" over all the tra- 
jectories starting from a fixed ipi and we will not integrate over this initial configuration. 
In the A, c, c we have put as indices (ft, ir to indicate either the first or the second set 
of indices "a" as in postulate II of section 2. To the currents we have coupled only the 
"configuration variables" in the representation (2.25) where <f>, X w , cP, c w are multiplicative 
operators. We did that because, as it is clear from (2.25), it will be straightforward to pass 
to the super-field formulation. The peculiar combination of "±" and "i" in the coupling 
of current and fields (3.3) is also related to the fact that eventually we also want to build 
a supercurrent to couple it to the superfield. 

As it is usually done, let us divide S into a free part: 



■ = j f dt J d 3 x X a <p a - A^tt - X n (V 2 - m 2 ) 



(3.5) 



+ic a c a - ic^c? - ic-x (V 2 - m 2 ) S 



and an interaction one: 



Sy = d x 



9 

3!^ 



(3.6) 



where we indicate with J d A x the expression f* f dt J d?x. To develop the perturbation 
theory, as it is usually done in quantum field theory, in the Sy of Z we shall replace the 
fields with the derivative operators with respect to the associated currents as follows: 



so we get: 



15 15 



1_S_ 

i 5J„ 



S 



j A„ 



1 5 



i 5J\ 



= / d 4 x 



5Jr 



9 $ 



5 3 ig 



8J, 



5 5 2 



3! 5J X „ SJl 2! SJs. 5 J 2 6J„ 



'<P 



i SJrf, i SJ^ 5J c <p SJ 5 
It is then easy to write JE^JJ^,, J^, J c <t>, JgJ in Eq. (3.3) as follows: 

"1 S 1 S 5 5 



(3.7) 



(3.8) 



•>, J\ n , J c 4>,JcJ = exp<iS v 



Ctt 



where Z$ is the free generating functional built out of the Sq. It is easy to prove that 

r(0) 



(3.9) 



this part Z^/ can be factorized into a "Bosonic" (B) part for the variables <p and A and a 

(3.10) 



"Fermionic "one (F) for the Grassmannian variables c, c: 

^Ipi [J^t J\k j Jc$ > Jo*] = Z {b) ipi \J<t>-> J>™\^\f) > Jc 
where 



Z (B) yS J <t» J *J 

= J 9'ipSiX exp j« jf 1 dt J d 3 x X a <j) a - A^tt - A,, (V 2 - m 2 ) </> + J^cp + J X7r X w j 



(3.11) 
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and 



Z (F) [^c*> J?*] 



= I @c@cexp{ I * dt I d 6 x 
hi 



ic a c a — ic^c 77 — ic w (V — m ) cr — iJ^c^ — ic^J-, 



(3.12) 

All the fields enter at most in a quadratic form. It is easy to integrate them out in Eq. 
(3.11) as it has been done in ref. [17]. The result is: 



A J = exp 



J d 4 xJ ( f,(x)(j)o(x) + i J d A xd A x' J^(x)Gr(x - x') J\„(x') 

where 4>o(x) is a solution to the equation 

(□ + m 2 ) 4> = 
and Gr is the retarded propagator satisfying 

(□ + m 2 ) G R = -5{x). 

Its expression is: 



Gr(x) = 



dp 



-ip.x 



(2tt) 4 (p° + ie) 2 - v 2 ~ m 2 ' 



(3.13) 



(3.14) 



(3.15) 



(3.16) 



Also for the we can easily integrate out the fields (for details see Ref. [30])and obtain: 



^^[J^.JgJ = exp 



d A xd A x' J c 4,(x)Gr(x — x')Jc n (x') 



(3.17) 



We have now all the tools to obtain all the propagators. 

Let us start at the zero order in perturbation theory and we will draw next to them 
the diagrams: 



<<K*)) C 



1 _s_ 

i 5J(i 

1 5 

i <5J^ <p. 



Z !pi [ J 4> > J K > J c<P > J cJ\ J^J^ , J + =0 



5 Z b\\ J ^ J A^]|j ,J A7r =O 



(f) (x) 



.1 



{(f)(x 2 )(f)(x 1 )) = 4>o(x 2 )4>o(x 1 ) = 



Xl 



X2 



ij SJ Xn (x 2 ) SJ^xi) 
= -iG R {xi - x 2 ) = G XlT <h = 



Z $[ J 4>i J \*, J c!t>-> J c*\\ 



(3.18) 



(3.19) 



(3.20) 



Xl 



X2 



In Eq. (3.20) we borrow from ref. [7]. the notation of full-dashed propagator . In Eq. 
(3.19) we have drawn a diagram made of two pieces not linked to each other, because they 
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are actually the product of two separated fields and not a propagator. We will have them 
soldered to each other once we use the temperature to average over the initial conditions. 
Analogously to Eq. (3.20) we can calculate 



{(f)(x 2 )Xn(x 1 ))° = -iG R (x 2 - xi) = G^ = •- 



x\ 



• 

X2 



(3.21) 



The correlations among Grassmannian variables give: 



(^(x 2 )^(x 1 )) c 



5Jc n {x2) SJ c4 ,(xi 



zP [J c0 , J c J|j j^=o 



(3.22) 



and 



{c < t>{x 2 )Cn(x 1 )Y 



Gr(xi - x 2 ) = Gc nC <t, = 

X\ x 2 



[J c <t> > J c n ]\ J , J E =0 



5J C 4,{X 2 ) 



(3.23) 



= -G R (x 2 - X]) = G c4 ,, = 



xi 



x 2 



Above we have adopted the convention of putting an arrow which points from c to c (see 
Ref. [30]). 

As both the propagators, of the (jxj) and cc, are related to the Gr we will sometimes 
put an index "(F)" or "(B)" to indicate if it comes from the Bosonic fields or "Fermionic" 
(Grassmannian) ones. The fact that the two propagators are equal (modulo i) to each other 
is due to the supersymmetry present in this formalism [11]. We will see other manifestation 
of it later on in several cancelations among diagrams. 

Let us now derive the rule for the vertices. Expanding e %Sv in Eq. (3.9) to the first 
order in g, we get : 



iS- 



v 



15 15 



5 



i SJs i dJs 5J r <t> SJ 5 



(fx 



1 6 5 3 



1 5 5 2 



216^6 J* 8 J* 



(3.24) 



Keeping account of the symmetry factors we get the following rules for the vertices : 




(3.25) 



where the continuos line refer to the field 0, the dashed one to the field A and the dotted 
one to c or c. 

Let us calculate the first order correction to the expectation value of the field (4>o(x)), 
i.e. (4>q(x)) (the super-index (1) on (</>} indicates the first order correction). Using the 
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vertices and the diagrams we get (where y is integrated over) 




The analytic expression of the first diagram above is 




- J d 4 y[^(y)G A ^(O)0 o (x)] (3.27) 
while the expression for the third one is 



x y 



= 9 J A [0o(y)<WO)M*)] . (3-28) 



x y 

Let us remember from (3.21) and (3.23) that 

G c ^ = -G R = -iG K4> (3.29) 
so (3.27) and (3.28) cancel each other and the only correction left in (3.26) is 

« ) = » < ( 3 - 3 °) 

whose analytic expression is 





ig J d 4 y[4>Uy)G4,xAx-y)] . (3.31) 



The reader could, at this point, claim that the loops of the retarded propagators 
appearing in Eq. (3.26) are both zero, because of the 9(t) which appears in all retarded 
propagators. For a loop to be zero we have anyhow to assume that 9(0) = 0. This 
regularization of the 6 corresponds to a prepoint choice in a discretized form of t [21]. The 
prepoint choice has some problems once it is used in the path-integral. It has in fact been 
proved in Ref. [23] that only if one uses the midpoint choice then the Feynman rules which 
can be read off the continuum version are the same as the one of the discretized version. 
So, for this reason, we have preferred the midpoint discretization which gives 6(0) = 1/2 
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and this implies that the retarded loops are not zero. Let us remember also that it is with 
the midpoint rule that the determinant is not a constant and the Grassmannian variables 
are needed. So, our formalism has an overall logical coherence. 

Going now back to our diagrams the reader may be puzzled by diagrams like (3.30) 
where there are just three lines at the end going nowhere or like the single line in (3.18) 
also going nowhere. This is due to the fact that our generating functional Z lfi depends on 
the initial configuration which is not averaged over. These kind of diagrams, but without 
Grassmannian variables, were first developed in [17] and are somehow similar to those de- 
veloped for fluid dynamics in [31]. Perturbation theory and something similar to Feynman 
diagrams where also developed long ago in [32]. We have not tried to compare our formal- 
ism with this other one. For sure our formalism is different from these other ones because 
it contains extra ingredients like the Grassmaniann variables. Moreover, because of these 
extra ingredients, we will be able to develop the super-diagrams which will simplify the 
diagrammatics . 

We will continue now to the second order for the (4>). The result that can easily be 
obtained (where y\ and are integrated over) is the following: 




As before, it is easy to see that there are cancelations. The second diagram is cancelled by 
the eighth one, the third by the seventh, the fourth is cancelled by the sum of the sixth 
and ninth, the fifth one by the tenth. So the final result is 




(3.33) 
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notice that no loop is left. 

The various symmetry factors have being obtained by working out the analytical details 
from the Z lpi in (3.9) up to second order. 

Let us now go back to the two-point functions which, to zero order, was given in Eq. 
(3.19). After some tedious calculations with the second derivative of JS^JJ], we get: 




(3.34) 



The first and the fourth diagram cancel each other and so the final result is 



W*2M*i)>« 



1 

3! x\ 



x 2 + 3! x\ 



X2 



(3.35) 



Going to the second order, after some long but straightforward calculations, we obtain: 



<#x 2 Mxi)>g> = 





_ 14 _ 







There are the usual cancelations between various diagrams and the final result is: 



{4>{x 2 )4>{ Xl ))^ 




Note that also here no loop is left. This is the real sign that we are doing a classical 
perturbation theory. When we will introduce the temperature later on loops will appear, 
but they are due to temperature and not quantum effects. Actually, loops would appear 
any time we do an average over the initial conditions so we can not strictly say that no 
loops is a sign of classicality, we can say that no loop is the sign of a generating functional 
which has no average over the initial configurations. 
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Collecting now all diagrams developed before for the two-point functions ,up to second 
order, we get: 



{<j){x 2 )^{x 1 )) ipi 



Xl 



1 

X2 o! Xl 




-• + 

X2 



1 1 

+ 2! 3! x*T 



(3!) 2 X! 




1 1 

~? 2 + 2! 3! x5~ 




^2 



+ 



X2 



(3.38) 



The rules for the symmetry factors are the following: 

1) put a factor of 1/n! where n are the number of "free" legs, i.e. those not starting or 
ending at x\ or X2] 

2) add a factor of 1/2! for any exchange of propagators which would leave the diagram 
invariant. 

Before concluding this section we would like to answer to a question that for sure many 
readers may have. The question is why all the diagrams in this section are disconnected or 
have legs going nowhere. The reason is not related to the fact that we may have used the 
wrong generating functional, but to the presence of the diagram (3.18). To convince the 
reader, let us give the analytical expression for some of the most "strange" diagrams. For 
example, the 10th diagram in (3.36) where there are both disconnected pieces and lines 
going nowhere: 



X2 



xi yi 



/ 



Ai d 4 y2Myi) 2 Mz2) 2 G { p(x 1 



X2) 



G\ \X2 - Xl) (/> (xi) (f)o(x2) 

We see that we have at xi and x 2 the fields (po(xi) and ^0(^2) and their diagram is the 
one of (3.18). In yi and y2 we also have two lines going nowhere and, according to (3.18) 
there must be two fields 4>o(y). The rest, between yi and 2/2, are two "Fermionic" retarded 
propagators. 

The reason we can not connect xi with yi or 1)2 with X2 with a continuous line is because 
we do not have this kind of propagator in the formalism with fixed initial condition. The 
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only propagators we have are those in (3.20), (3.21), (3.23) and (3.22), which are the dash- 
full propagator or the Fermionic dot propagators. For the <f>(f> we have only the disconnected 
propagators (3.19). The reader may think that by coupling in some way various dash-full 
propagators one could get a (f><f> propagator. If one tries, it is easy to convince himself that 
this is not possible. 

In this chapter we have studied the " average" of a single field <p or correlations of two 
fields 4>4>. It is possible to do the same for the fields and the c^, c w ones. Instead of 
doing them explicitly we will work out the super-field perturbation theory for this Z lf>i [J] 
in appendix A. We can then project out the components of the super-fields and get all the 
correlations mentioned above. We advise the reader not to jump immediately to appendix 
A, but to wait till he has read sec. 4. 



M*,t) = I 7^sMP,t)e^ (4.2) 



4. Perturbation Theory With Temperature. 

Let us now suppose that, instead of working with some fixed initial configuration, as in 
section 3, we do a thermal average over the initial configuration. Let us use the following 
notation: x for the 3-dim vector and x for the 4-dimensional one. Let us remember that 
in section 3 we choose a solution 4>o(x,t) of the Klein-Gordon equation: 

(df- V 2 + m 2 )0 o (x,t) = 0. (4.1) 

Its Fourier transform <fio (p, t) is : 

d 3 p 

and satisfies the equation: 

{d 2 t - E}) MP,t) = 0, (4.3) 

whose solution has the form 

MP, t) = UP) cos [E f (t - U)] + ^ sin [E p it - t t )] (4.4) 

with ti, 4>i, TTi the initial time and field configurations. So the final expression for the field 
is: 

M*, *) = / (0» [h® cos " *<)] + ^ sin W ~ *<)] • ( 4 - 5 ) 

The thermal-averaged correlation- functions that we want to evaluate are defined, for ex- 
ample for a two-point function , as follows: 

/^(f)^(x)(#) WE(ilj) e^*^ 



/ S>(j)i{x)SiT:i{x)e-P H ^i) 



(4.6) 



where (f>i(x), 7Ti(x) are the Fourier transform of the <j>i(p) and n(p) appearing in Eq. (4.4). 

On the correlation {(pfyp we have not indicated the argument of the fields. We did 
that in order to simplify the notation for the moment. 
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The analog of the generating functional (3.9) will now be: 



(4.7) 



In order to do the integration over the initial configuration in (4.7), let us first do the 
anti- Fourier transform of the 4>(p) and tt(p): 



So we can write <po(x,t) as: 

4>o(x, t) = J d 3 af(f)i(x')a ( j)(x — x' ,t) + J d 3 x'iri(x')a n (x — a?, t) 

a^x - x' , t) = J dPxe'^ 3 -^ cos [Ep(t - U)] 

J Ep 



where 



(4.8) 
(4.9) 

(4.10) 

(4.11) 
(4.12) 



Next, let us extract from (4.7) the part which depends only on the initial field config- 
urations which is: 

J ^ i (f)^7r i (x)e-^ ^' ,ri > + J" d4a!j *^* W 
where we have switched off [10] the interaction in H so that Hq is just: 



#o(<&,tt;) = J d 4 



n (Vc^) 2 m 2 tf 
2 2 2 



(4.14) 



Inserting (4.10), (4.11) and (4.12) in (4.13) it is easy to perform the integration over the 
initial condition (see Ref. [10] for details or appendix A of our paper) and the result is: 



Zp[J$] = exp 



X - J d^xd^x'J^Afsix-x')^') 



where: 



It is easy to prove that 



Ap(x-x') = Ap(x'-x). 



(4.15) 

(4.16) 
(4.17) 
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Going now back to the full expression (4.7) of the generating functional, the full Bosonic 
part is: 

z {B)p[J^] = ~\ j d 4 xd A x'J 4> (x)A /3 (x-x')J ( f > (x') 

+i J d 4 xd 4 x' J ( f,(x)G R B \x — x')J\ iv (x') 
Let us calculate the thermal two-point function 

5 



(4.18) 



(4>(x 2 )4>(x 1 )}i 3 



5J < f > (x2) 5J<j,{xi 



(4.19) 

x 1 



We can see that this correlation is not anymore the product of the two fields at x 2 and 
like in (3.19), but it is a function Ap(x 2 — x\) which links x 2 with x\. We will use the full 
line to indicate the Feynman diagram associated to Ap{x 2 — x\) 



Ap(x 2 - xi) 



X\ 



x 2 



(4.20) 



Let us now see how to get the first order correction to the two-point function. Let us, for 
example, look at the first diagram in (3.35), and let us take its thermal average indicated 
by ( )p: 



3! 



Xi 




x 2 



= Jj J dyG R (x 1 -y)((f>o(y) 3 (j)o(x2)) t 

= § / d y°R^-y) \ 



4 5J^{yf5J^x 2 ) 
- j dyG R {xi-y) 3Ap(y,y) Ap(y — x 2 ) 



J<4=0 



1 

2! xV 



x 2 



(4.21) 



Let us notice that we do not have anymore disconnected diagrams. The tree legs in y going 
nowhere and the one in x 2 get soldered to each other in all possible combination producing 
the loops in (4.21) in y and the propagator between y and x 2 . 
The same can be done for the second diagram in (3.35) 




(4.22) 



The same for the second order contribution in (3.37). After some long but straightforward 
calculations, we get from the first diagram in (3.37): 
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-- - 

2! 3! \ Xl 



X2 



2! 



xi 



X2 



(4.23) 



(2!) 2 X! 



+ 



x 2 2! xi 




For the second diagram in (3.37) we obtain: 



ii - 

2! 3! \ xi 



^2 



2! 



+ 



xi 



X2 



(4.24) 



(2!)2 Xl 




— • 

^2 



and the third term in (3.37) gives: 



(3!) 2 \ Xl 



X2 



(2!)3 xi 



-• + 

^2 



3! xi 



Summing up all the terms up to 2nd order we get: 




X2 



(4.25) 



(</>(xi)<t>(x 2 ))p 



Xl 



x 2 2! xi 



x 2 2! xi 



-• + 

X2 
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2! x\ %2 2! x\ %2 




(4.26) 

The reader may wonder that there could be some extra diagram coming from the Fermionic 
part of the generating functional that we have not considered. This is not so because in this 
calculation we started from the result we had got for the analog diagrams without thermal 
averages (3.38) and for those diagrams we had already taken account of the Fermionic 
diagrams cancelation. For example, just using topological considerations we could have for 
the two-point functions also the following diagrams besides the (4.26) 




but all these would be canceled by the diagrams with Fermionic loops like: 
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while no one of those in (4.26) could be canceled by diagrams with Fermionic loops. 

The diagrams we get by averaging over the initial conditions with the Boltzmann 
weight are basically based on the following propagators and vertices 



2tt 



5(p 2 — m 2 ) 



Xl 



P 



2 — m 2 + iep° xi 



iGf(-p) 



Xl 



-ctP(p) = -Gg>{-p) = 

X\ 

ig 



X2 



• 

X2 



X2 



Xl 



X2 



■ m 

X2 




(4.28) 
(4.29) 
(4.30) 
(4.31) 
(4.32) 

(4.33) 
(4.34) 



5. Perturbation Theory With Superfields. 

If we add to the of (4.18) the Grassmannian piece (3.12) (which is not affected by the 
thermal average) we get the full generating functional which, for the free theory, is : 



■\ J <fxd A x l J 4> {x)/\p{x-x')J l p{x') + 
+i J d 4 xd 4 x / J (/) (x)G { p(x-x / )J Xn (x')+ 
-i J d A xd A x' J c4> {x)G ( P {x - x')J- c ^x' 



(5.1) 



Let us now move quickly to the superfield formalism. The super-field analog of (2.25) is 
the following one: 

$ = + Qc^ + 0c* + i96K (5.2) 



and in order to have the current-field coupling that we have in (3.11) and (3.12), we need 
to introduce a super-current of the form: 



-j^ + ej-^ + oj^-mj^. 



(5.3) 



The term: 



cxp 



j d 4 xd9d9j (f> (x,9,9)$(x,9,9) 



(5.4) 



will then give the four couplings of currents and fields present in (3.11) and (3.12). If we 
want to get not the free generating functional but the full one, its formal expression will 
be 

5 



I 



dzV 



Z[3$] = exp 
where V is the potential of the g^ 4 /4! theory, so 

vt- 5 ^ o sl 



SMz) 



(5.5) 



SMz) J 4!5J4(z) 




g J d 4 yd9d9 = g j dz, 



(5.6) 



where the collective variable z = (x, 9, 9) has been introduced. The reason we can use this 
potential and not the <Sy of Eq. (3.10) is because, as we proved in Ref. [19], the Sy or any 
<S can be obtained via the usual S with <j) replaced by $ as in formula (2.18) and (2.19). 
A super-field formalism and the associated Feynman diagrams had been developed for the 
Langevin stochastic equation in ref. [36] . This has inspired us to do an analogous formalism 
for deterministic systems which evolve via Hamilton's equation. 

Going now back to the coupling between supercurrents and superfields in (5.4), it is 
straightforward to prove that: 



J (f) = i J d9d9J4x,9,9) 

J x „ = - J d9d999J (f) (x,9,9) 
J c * = - J d9d93 lj) {x,9,9)9 



J 5n = / d9d99J4x,9,9). 



(5.7) 

(5.8) 
(5.9) 
(5.10) 



Using the rules (5.7), (5.8), (5.9), (5.10) and the super-field we can write the free generating 
functional (5.1) as 



Z2[J ] =exp 



- / dzdz'j^z)G{z,z')3^z') 



(5.11) 



where z = (x, 9, 9) and G(z, z') is defined as: 



G(z,z') = Ap(x-x') + g{z,z') 



(5.12) 
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with 

G(z, z) = G^\x - x')6'6' + G ( k\x' - x)60 + G% \x - x')60' + G { p{x' - x)9'e. (5.13) 

The G(z, z') and Q(z,z') have several nice properties that are the following (with m > 1) : 
1. 

G(«,«') = G{z',z). (5.14) 



2. 



G(z,z) = Ap(x-x). (5.15) 
g n (z,z') = <5 n ig n (z, 2! / )[l- 5(2- 1 (5.16) 



4. 



G m (z,z') = A^ , (x-i')+mA™- 1 (a;-i')^^'). (5-17) 

While 1) and 2) are trivial to prove from the symmetry properties of and Q, 4) 
can be proved as follows: 



G m (z,z')= f m ) ^' n (x-x f )G n (z,z') 

0<n<m \ / 



Af{x-x')+ I - I Aj-\x -x')g(z,z') (5 " 18) 







(5.19) 



= A™(x - x') + m AJ- Y {x - x')G(z, z'). 

In the last step we have used Q n (z, z') = 0(n > 2) which is a consequence of the property 
number 3). Let us now prove property number 3). We have, for n > 1: 

g n (z, z') = G n (z, z>) {5{z - z') + [1 - S(z - z')]} 

= g n {z, z)5(z - z') + g n {z, z')[l - 5{z - z')} 
= [S nl + (l-S nl )]g n (z,z')[l-6(z-z')] 
= S nl g n (z,z')[l-d(z-z')} 

+(l-Sni)g n (z,z')[l-5(z-z>)} 
= 5 nl g n {z,z')[\-5{z-z>)]. 

The property no. 3) is at the root of the fact that loops made of dash-full line cancel against 
Grassmannian loops. So it is at the root of "classicality" , i.e. that without temperature 
we would not have any loops. 

Going back to the expressions (5.5) and (5.11) we have seen that the Feynman diagrams 
are the same as those of a g^/Al theory, but with the field replaced by super-field, and the 
propagator replaced by the super G(z, z'). We call the reader's attention that, in contrast 
to what happens in standard quantum field theory, where, due to translational invariance 
we have Gp{z, z') = Gp(z — z'), where Gp stays for the Feynman propagator, here we have 
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G(z, z') / G(z — z'). That is why we use the notation G(z, z') in the perturbative analytic 
expressions. 

Let us start deriving the two-point function: 



G f) {z l ,z 2 ) = {<S>{z l )<S>{z 2 ))p = {-i) 



■Z[J«] 



JU=0- 



(5.20) 



Doing the analog of standard gcp 4 /A\ QFT perturbation theory, we have that the first order 
correction to the two-point function is 



X - f dzG(z 1 ,z)G(z,z)G(z,z 2 ) = \ 



2 z x 



Z2 



(5.21) 



which is the super-analog of the QFT diagram with the same symmetry factor. The second 
order correction, following the analogy with QFT, which has the three diagrams: 



3! 




-• + 



-• + 




is made of the following three super-diagrams: 




i J dzdz f G( Zl ,z) [G(z,z')] 3 G(z',z 2 ) (5.22) 




dz dz'G{zi,z) G(z, z) G(z, z')G{z , z) G{z , z 2 ) 



(5.23) 



=] f dzdz'G(z 1 ,z)G(z,z')G(z',z')G(z',z)G(z,z 2 ). (5.24) 

z 1 z z 2 4 J 

Let us now check if from the super-diagram (5.21) we can derive, for example, the 
correlation at the first order which were the sum of (4.21) and (4.22). The manner to do that 
is to extract from the external super-fields of Eq. (5.21) the components <p a . Remembering 
the formula of the super-field (5.2), this can be done via the following integration 
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J MxdBx (Mi) d6 2 d6 2 (M2) 



2! Zl 



z 2 



(5.25) 



The variables inside the brackets and (0 2 6 2 ) basically act as projectors from the 

super-field <1> to the field <fi. 

Let us now work out (5.25) in details using (5.21) and the analytic expression of the 
super-propagators presented on the l.h.s. of (5.13). Below we will go pedantically through 
the details of the calculations. We do that in order to get the reader familiar with the 
formalism. 

The explicit expression of (5.25) turns out to be 




J d0 1 d9 1 (Mi) d0 2 $2 (M 2 ) 
= \ J dO-tdh (Mi) dd 2 d0 2 (M2) J d 4 xd9d6 

x Ap(xi -x) + G { *\x - xi)Mi + G { *\ Xl - x)6e + G ( p(x - xi)00i + G { p{x 1 - x)6 l 6 
x [Ap(x-x) + g(z,z)] 

Ap(x - x 2 ) + G^\x - x 2 )9 2 2 + G^\x 2 - x)60 + G { £\x - x 2 )B~B 2 + G ( p{x 2 - x)9 2 e 

Ap(x - x') [Af,(x - x 2 ) + Gp\x 2 - x)0e 
&p(x - x') 



= \ J d 4 xd6d9 Ap( Xl -x) + G { *\ Xl - x)66 
= \ J d 4 xd6d6 A /3 (xi -x) + G { n\ Xl - x)06 
Ap(x - x 2 ) + G ( P (x 2 - x)0e 



(5.26) 



1 

2! x*T 



x 2 2! x\ 



x 2 



and these last diagrams are exactly the sum of (4.21) and (4.22) which were the first order 
correction to the (4>4>) propagator. So the lesson we learn from here is the following: it is 
enough to project the external super-legs of the super-diagrams on the fields we want, in 
order to get also the correct internal part of the super-diagrams. 

Let us now check the second order. That means let us check if the three super-diagrams 
in Eqs. (5.22), (5.23) and (5.24), once the external legs are projected on (j>, reproduce the 
eight second order diagrams contained in (4.26). 

Again, in order to get the reader familiar with this formalism, let us calculate explicitly 
here the first of the three super-diagrams in (5.22), while the calculation of the other two 
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will be confined in appendix B 




(5.27) 



JdB 1 S 1 de 2 S 2 (Mi) (M 2 ) 

= J dd x dQ x d9 2 dd 2 (Mi) (M2) J 



3! ^ z 2 2 

d6 1 d6 1 d9 2 d6 2 (Mi) (6 2 6 2 ) J d 4 xd9d9 J d A x'd9'd9' 

A p ( Xl -x) + G^\ Xl - x)99 + G { x\x - xi)Mi + G { p{xi - x)9 1 + G { p{x - x 1 )99 1 x 
[Ap{x-x')+g(z,z')] 3 x 

Ap(x' - x 2 ) + GSf \x' - x 2 )9 2 9 2 + G^\x 2 - x')9'9' + G% \x' - x 2 )9'9 2 + G^\x 2 - x')9 2 9' 
= j d9 1 d9 1 d9 2 d9 2 (Mi) (9 2 9 2 ) J d 4 xd9d9 J d 4 x'd9'd9' x 

A^xi - x) + G% - x)00 + Gjf^x - xi)Mi + G^f } (xi - x)M + G^f } (x - xi)00~i 
A^(x - x 1 ) + 3A 2 (x - x')G(z, z' 

= -J d 4 xd 4 x'A l3 (x 1 - x)A§(x - x')G { *\x' - x)Gf\x 2 - x' 
+-J d 4 xd A x'G ( p(x 1 - x)A|(x - x')G { *\x - x')Ap(x' - x 2 ) 
+-J d 4 xd 4 x'G ( ^\x l - x)A^(x - x')G^\x 2 - x') = 



+ 




1 

+ 7^7 



x 2 2! xi 




+ (5.28) 



3! xi 



^2 



These are exactly the last three diagrams of (4.26) and exactly with the same coefficients. 
Analogously we can prove (see appendix B) that 



J d9 1 d9 1 d9 2 d9 2 (Mi) (M2 



(2!) 2 zi 





z 2 



(5.29) 



(2!) 2 X! 



+ 



x 2 (2!)2 Xl 



+ 



x 2 
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+ 



(2!) 2 *i 



X2 



which are the 6th, 7th, 8th diagrams presented in (4.26). 
Along the same lines we can easily get (see appendix B) 



J d6 1 de 1 de 2 d6 2 (Mi) (0202) 



(2!)* 



1 

2! x*T 




(5.30) 



X 2 A Xl X 2 

Summing up (5.29), (5.30) and (5.27) we get all the second order diagrams contained in 
Eq. (4.26). These last ones were eight in numbers and we got the from just the three 
super-diagrams in Eqs. (5.22), (5.23) and (5.24). This decrease in the number of diagrams 
is only one of the virtue of the super-diagram technique. A second one is that we only 
have to perform integrations over the Grassmaniann variables (which are easy to do) to 
get all the old diagrams of L without bothering with their complicated symmetry factors, 
without bothering with soldering in the correct manner the legs and vertices of £ and so 
on. Above all, the main advantage of the super-diagrams is that they can be derived from 
the analog quantum theory whose Lagrangian L is in general simpler than C. So, once we 
have the quantum theory and its associated Feynman diagrams and symmetry factors, we 
are already more than halfway also with the perturbation theory of the associated classical 
theory. We have just to substitute the fields with the super-fields and the propagators with 
the super-propagators and perform, as we said above, some very simple Grassmaniann 
integrations. Moreover, the super-diagram automatically perform for us the cancelation 
between Fermion and dash-full loops. 

Last, but not least, by changing the projectors on the external legs, we can get not 
only {(pip) correlations, but also {\ip), the {cc), the (cc), the (cc) and the (AA) correlations. 
So just three second order super-diagrams will produce tens of standard diagrams. 

We will start now with the {(j)X w ) correlations. Looking at the super-field expression in 
(5.2), we realize immediately that, differently than the (</></>) correlation, we have to project 
out only the 4>{xi) field because the X n {x 2 ) is already equipped with its own 0202 an d so it 
naturally makes in appearance once we integrate over the final points. From the vertices 
and propagators that we can build from the formalism in components, it is easy to see that 
{4>\ n ) correlation will have the zero-order component of the form: 



{4>{xi)K{x 2 ))o 



Xl 



• 

X2 



(5.31) 
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While the first order correction is: 

(ip(x 1 )X n (x2))i = [d0 1 d§ 1 d6 2 $2 {-ididi) 



1 

2! xT 




X2 



1 

2! 



^2 



(5.32) 



Here we have used the projector (— i6i6\\ in order to project 4>{x\) out of the super-field. 
The second order corrections have the form: 



(<j){x l )K{x 2 ))2 = Jd0 1 dB 1 d0 2 S 2 {-ihOi) 




+ (5.33) 



(2!)2 Zl 



By long but straightforward calculations similar to those performed in appendix B, the 
result that we get for the diagrams in (5.33) is the following one: 



3! xi 



+ 



2! 



X\ 




+ 



X 2 (2!) 2 X! 




(5.34) 



We checked and found out that the diagrams on the r.h.s. of (5.34) are exactly those 
that we would have obtained using the component formalism with the rules given in Eqs. 
(4.28), (4.29), (4.30), (4.31), (4. 32), (4.33), (4.34) and taking into account of the cancelations 
between dash-full propagators and Fermionic ones. 

From now on we will trust the results given by super-fields once they are properly 
projected on the components fields we are interested in. We will trust them because they 
seem to give the same results as if we have done the calculations in components using the 
rules in Eqs. (4.28) (4.29), (4.30), (4.31), (4.32), (4.33), (4.34). 



- 29 - 



Let us now calculate the Grassmaniann variables correlations. That means we consider 
them as physical fields and not like "ghosts", as it is done in gauge theories, where they 
contribute only to internal loops and are never physical (i.e. external ) fields. The reason 
we consider them as physical fields is because, as we said before, they are the Jacobi fields 
and they enter in the calculations of the Lyapunov exponents [24]. 

Let us start from the correlation (c^ (x\)o K {x2)) ■ The zero order is given in Eq. (4.30). 
The first order, via super-fields , will be obtained from the super-diagram of Eq. (5.21) by 
projecting out the c^(xi) and c w (x2)- This is achieved with the projector #2#i- The first one 
extracts from the super-field the field c 7r (x2), while the second one extracts the c^{x\). So 



1 



what we have to do is 

2! z\ z Z2 

Performing the long but straightforward Grassmannian integration above , we get: 



J de 1 de 1 de 2 de 2 (e 2 9i) 



i 

2! 



Xl 



X2 



(5.36) 



For the second order we just have to perform the following calculations: 




We leave this calculation to the reader and we let him compare it with the calculations 
in components which we can get from the rules stated in Eqs. (4.28), (4.29), (4.30), (4.32), 
(4.33). 

Next, let us calculate 

{c't > {x 1 )^{x2))p. (5.38) 
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The projector in this case is (#16*2): 

J d9 1 d9 1 d9 2 d9 2 {hh) ^{zt)^{z2))p = {c^{x 1 )c ( l > {x 2 ))p. 



(5.39) 



By an argument presented in appendix D we can prove ,to all orders in perturbation theory, 
that: 

{c*( Xl )c*{x 2 )) ri = 0. (5.40) 
With similar arguments we can also prove that 



{x{]c 4> {x 2 ))p = {(t>{xi)o K {x 2 )) ^ = (A(xi)c^(x 2 ))/3 = (A(xi)c 7r (x 2 )) ( g = 



(5.41) 



Let us just work out the first one of the relations in (5.41). The projector needed to 
derive things from the super-fields is {9 2 9\9i), so: 



d9 1 Side 2 dS 2 (0 2 Mi) ($(zi)$(-z 2 )>/j 

J d9 1 d9 1 d9 2 d9 2 {O2O1O1) J dzdz'G(z 1 ,z)¥^(z,z')C(z',z 2 ) 



I 



(Wi) J , 



d9 1 d9 1 d9 2 d9 2 (l),(h(h) / d 4 xd9d9 I d A x'd9'd9' 
(B) (n , _ ^ ^ a + q(F) ^ _ ^ p $ n (F) 



Ap( Xl -x) + G { ^\ Xl -x)00 



+g { k\x - x 1 )e 1 e x + G { p{ Xl - x)e 1 e + G { p{ x - x 1 )ee 1 fW( Z) z ') 



(5.42) 
(5.43) 



x Ap(x' - x 2 ) + G^>(x' - x 2 )9 2 9 2 + G K °\x 2 - x')9'9' 
+G { p{x' - x 2 )9'9 2 + G% \x 2 - x')9 2 9' 
= A + B. 

where F^ ra ^ is defined in appendix C, while A and IB are: 

A = J dzdz'e'Ap( Xl - x)¥ < - n \z,z')G { p(x 2 - x') 
and 

B = / dzdz'm'G { p( Xl - x)¥ {n \z,z')G ( p{x 2 - x'). (5.44) 

In appendix C we proved that the Grassmannian coefficients of F^ n ^ can only be one of the 
following six forms: 1) 999' 9' , 2) 90, 3) 9'9', 4) 99', 5) 99' and 6) D. Inserting any of them 
in A or B we get that A = = B in all cases and this proves the first of the Eq. (5.41). The 
proof of the other relations in Eq. (5.41) are analogous and we leave them to the reader. 
Let us now prove that, to all orders in perturbation theory, we have: 



(A 7r (xi)A 7r (x 2 )) / 3 = 0, 



(5.45) 



To simplify the calculations let us use the notation d/x = d0d0. In order to project out the 
A field from the super-field we need, as projector, the operator —I, so 



/ 



dM0ld02d02(-n)($(*l)$(z2)}/J = (K(xi)\Ax2))p. 



(5.46) 
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Let us now consider (5.46) at n-th order in perturbation theory 

J d0 1 de 1 de 2 d0 2 (-\[)(^(z 1 )^(z 2 )}P 

= - j d9 1 d9 l d9 2 d9 2 j dzdz'G{z 1 ,z)¥ { - n) {z,z')G{z',z 2 ) 



- I d9 1 d9 l d9 2 d9 2 I d 4 xd9d9 I d 4 x'd9'd9' 



Ap(xi -x) + G ( x\xi -x)99 



+G ( k ] (x -x 1 )e 1 0i + G ( P (xi -x)9 1 9 + G ( P (x - xi)00i] (z, z') 
x Ap(x' - x 2 ) + GjfV - x 2 )9 2 9 2 + G%\x 2 - x')9'9' 
+G ( p(x' - x 2 )9'9 2 + G { p(x 2 - x')9 2 9' 

= - [ dzdz' [ d9 1 d9 1 d9 2 d9 2 \g^\x - x^O^] ¥^ n \z,z') - x 2 )9 2 9 



(5.47) 



= - j d 4 xd 4 x'G { p(x - xi)GSfV - x 2 ) j dndn'¥^ n \z,z') 
In appendix C we will prove several identities which will lead to 

dfidfi'¥ {n) (z,z') = 0. 

Because of this, from (5.47) we obtain, to all orders, that : 

(A 7r (xi)A 7r (x 2 )) ( g = 0, 



(5.48) 



(5.49) 



The correlation above is zero not only because of the identies that we proved but also 
because we somehow are using here the closed-time-path (CTP) formalism for classical 
thermal field theory [30]. We will not expand on this here but advice the interested reader 
to read ref.[30]. 

For the correlation {<j)(x 2 )<j)(x\)) p we have calculated things up to second order which 
is what has been done in the literature up to now. We will show in appendix E that with 
our super-field technique it is not difficult to reach the third order. 

6. Fluctuation-Dissipation Theorem. 

This is a well known theorem that holds for many systems like those obeying the Langevin 
equation and similar equations. It basically relates the two-point function (fluctuation) 
of the system to the manner the system responds to external perturbation (dissipation). 
It was proved to hold [16] also for deterministic systems whose initial distribution was 
the Boltzmann one, so it should hold also for the approach to classical field theory with 
temperature that we have explored in this paper. 

Several years ago it was showed [33] that the fluctuation-dissipation theorem (FDT) 
for Langevin equation could be proved non-perturbatively as a "Ward identity" of a hidden 
supersymmetry present in the Langevin equation. As that supersymmetry is very similar 
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to the one we have in the classical formalism presented here, also in the classical case the 
FDT can be derived as a Ward identity of the supersymmetry. In appendix F we will 
present that non-perturbative derivation generously provided to us by Martin Reuter. Of 
course, the FDT for deterministic systems can be proved in a much simpler way [16]than 
going through the Ward identities of susy, nevertheless, we like to look at it as a Ward 
identity because it actually relates different correlation functions like Ward identities do. 

For the moment , anyhow, let us leave aside its relations to Susy and Ward identities 
and let us start verifying its validity at the perturbative level, so that we can test the tools 
developed in this paper. Let us follow the notation of Ref. [16], that from now on we will 
indicate with the initial letters "DH" , while we will indicate ours by the acronym CPI. In 
the table below we write down the correspondence between the notation of DH and the 
one of the CPI. We urge the reader to read Ref. [16] in order to understand the symbols 
contained in the table below: 



[DH] 


[CPI] 


Dij 




Rij (t — t') 


-ie(t-t')( v a (t)\ b (t')) 


Cij{t-1!) 


(<P a (t)<p b (t>)) 






^ = s$~i 





The ip a indicates all the phase-space coordinates. According to Eq. (2.34) of DH the FDT 
is 

d_ 

di' 



(6.1) 



where we follow the convention adopted in DH: the bar over an index, like k, means they 
are summed over. 

Eq. (6.1) can be written in an alternative form as 



Q \ 1 

— Cij(t) = -R ik (t)D kj - -D ki R jk (-t). 
We can prove it in the following way. For t > 0, Eq. (6.1) becomes 



(6.2) 



Rii(t) = -P^C fk D- kj . 



(6.3) 



If we multiply the relation above by D we get: 



1 d 

-pRili^Dtj = --Q t C ik D ki D lj- 



(6.4) 



Next, using the fact that D k iD 



-5kji we obtain: 



8 C ij {t) = - p R ik {t)D- kj . 



dt 



(6.5) 
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For t < Eq. (6.1) gives 



d_ 

= m-t)^-c fk (t)D- kl 



(6.6) 



where we have used the fact that 



§-fi k (t) = ~C kj (-t). (6.7) 



Multiplying (6.7) by D and operating as before, we get for t < 



± Cij {t) = -±D ki (t)R jk (-t). (6.8) 



Collecting (6.5) and (6.8) for t ^ we get exactly (6.2). 

Making use of the table of comparison between DH and CPI, we can turn (6.2) into 
the following relation 

A^)/ (i 2 )) = ±0(ti - t 2 )(<p a (t 1 )(-i)\ d (t 2 ))u db - ±9(t2 - t 1 )a; da (/(t 2 )(-,)A (J (t 1 )). 

(6.9) 

This is the full fluctuation-dissipation theorem. 

Let us choose a = b = 1 in Eq. (6.9). We will then obtain for a field theory, where 
(p 1 = 4> and ip 2 = 7T, the relation : 

A^«( ti y (* 2 )) = i _L [(9 ( tl _ taX^A^ta)) - 0(t 2 - ti)Mt 2 )M*i)>] . (6.10) 

At order zero in perturbation theory, in terms of Feynman diagrams, Eq. (6.10) can be 
written as follows: 

q Ap(x 1 -x 2 ) i G ( ^\x 1 -x 2 ) i G { x\x 2 - xi) 

jrr • • = -3 • • - -s • • (6.11) 

Ot\ X\ X2 P X\ X 2 P X\ X2 

We will provide a proof of this in appendix G. 

At the first order in perturbation theory the diagrammatic form of the FDT (6.10) is: 



Id \/ Id W , 

— — — • V • • V • = (6.12) 

2 oti x\ x 2 2 ot\ x\ x 2 



It W 1 i 

• V • 



2 /3 x\ x 2 2 /? xi x 2 



and also this will be proved in appendix G. 
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The second order result is given by the following three diagrammatic relations. 




_L d_ 

(2!)2 0ti xV 



1 d 

-«+-TT-^— •- 



x 2 {2\fdh Xl 



-• + 

X2 



_L 9_ 



—m 

X2 



(6.14) 




2! ati 



+ 



1 d 



x 2 2! <9ii xi 



X2 



(6.15) 
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We feel that all these diagrammatic identities will turn out to be very very usefull in 
simplyfing long sum of diagrams which may appear in phenomenological applications of 
our formalism. 



7. Conclusions and Outlooks. 

In this paper we have developed the Feynman diagrams for a classical g - ^ field theory 
either with and without temperature. This topic is of high interest in the field of heavy 
ions scattering where classical field theory with temperature seems to be playing a central 
role. Besides this application, the Feynman diagrams for a classical theory can be useful 
in many fields from planetary motion, to fluid dynamics, to chaotic motion, etc. We have 
been able do develop the Feynman diagrams in a very easy way because a path integral for 
classical systems had been previously developed [11]. With respect to other approaches to 
perturbation theory for classical systems, we have been able to give the Feynman diagrams 
not only for the Bosonic variables of the system and for the response fields, but also for the 
Jacobi fields associated to the dynamics. Actually are these last fields (the Jacobi ones) 
which can provide us with indications on the cahotic behavior of the system. Our Feynman 
diagrams could for example be used to calculate perturbatively the Lyapunov exponents 
of a dynamical system. 

Besides developing this diagrammatics for all these fields and their interactions, we have 
showed that the many different classical diagrams of the various different fields mentioned 
above can be unified in few super-diagrams which have the same kind of vertices and kinetic 
term as the quantum one associated to ip a . 

We hope that this super-diagram formalism can be of some help not only in simplifying 
the notation, but also in understanding the subtle interplay between the quantum high 
temperature behavior and the classical one. Interplay which is very important in the 
heavy ions scattering field. We did not address this last issue here because this paper is 
only aimed at giving the formal diagrammatics necessary in this field. 

What we will do in a forthcoming paper [35] is to develop a formalism analog to the 
one of the "effective action" and try to give something like a " renormalization group" 
approach. All this is done in the spirit of provinding the formal tools which could later be 
used by physicists for more phenomenological applications. 
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A. A Fundamental Identity. 



In this appendix we will show the details on how to go from (4.13) to (4.15) If we insert 
(4.5) into (4.13), we get: 



I 



d 6 x 



^(V^) 2 , m 



2^2 



+ 



2 2 

xi / d 4 xJ^(x) / d 3 x4>i(x')a^,(x — x ,t) 



'I 



x j ^7rj(x)exp 

xi j d 4 xJ ( / ) (x) J d 3 x'ni(x)a 7T (x — x' ,t) 
Let us now define the following quantities: 

J^(x) = i J dtd 3 xJ ( f ) {x,t)a ( j > {x — x ,t), 

J%(x') = i j dtd 3 xJ^{x,t)a n {x-x',t) 
and let us rewrite the normalized Zp[J<j>] using (A. 2) and (A. 3): 

z p[ J <t>] 

(V<f>i) 2 | m 2 <f> 2 



where 



+ J d 3 f^(x)j;(x) 

X ^ / ^TTi(x)< 



exp <^ -/3 / drx 



+ / d 6 xiri(x)J%(x) 



x) exp < — / d x 



(V^ + m^2 



^[0] 7r = y^7r i (x)exp|-/3 y d : 



21 



Performing a partial integration for the kinetic piece in Eq. (A. 4) we get : 



exp {~2 d 3 xd 3 x'4>i(x) (36(x - x') ( 



■V 2 ' + m r 



(A.l) 



(A.2) 
(A.3) 



(A.4) 
(A.5) 

(A.6) 
(A.7) 



)] 4>0) + f d 3 x4>i(x)^(x) 



(A.8) 

As the weight is quadratic, the integration can be formally done for all (A.4) and the 
numerator of (A.4) turns out to be: 



exp 



(A.9) 



-37- 



where A(x, x') = f38(x — x'){— V 2 ' + to 2 ) and it obeys the equation: 
/ d 3 x 'A(x, x^A- 1 {x 1 , x") = 5(x - x") 
= f d 3 x" \l35(x - x')(-V 2 ' + m 2 )] A-\tf, x") = 5(x - x") 



so 



(A.10) 



(A.11) 



[/3(-V 2 + m 2 )A~ l ] (x, x") = 8(x - x"). 
Performing the Fourier transform of A" 1 we get: 

A^&x") = ^3 J cPpe-^'-^A^ip). 
It follows from (A. 10) that 

A = /3(^ + m 2 )' 

Using this and the expressions (A. 2), (A. 3), (4.11), (4.12) we obtain, for the argument in 
the exponent of (A. 9,) the following expression: 

J d 3 xd 3 x'^i J d 3 XidtiJ^(xi,ti) J ^^ e #i-(*i-2) cos [E Pl (h - U)] 



(A.12) 



(A.13) 



J (2vr 



d 3 P jp.tg' 



/ / d 3 x 2 dt 2 J <p (x 2 ,t 2 ) I 1?^ COB [Efrfa-U 



(2vr) s 



(A.14) 

where ti is the initial time appearing in (4.4). Performing, in the expression above, the 
integration over x and x' we get: 

- J d 4 x 1 d 4 x 2 J^x 1 )J4x 2 ) J -j^Le-V^-^A- 1 ® cos [E^h - U)] cos [E p {t 2 - t t )] . 

(A.15) 

So (A.4) becomes: 



x) exp 



(V^) 2 | rn 



-13 J d 3 x 

= exp |-i J d 3 xd 3 x'J^(x) | J jLLa- 1 ® cos [E^t - U)] 
x cos [E p (t' - U)] e ^'- £ )} J^(x')} 
Let us now do an analog set of manipulations for (A. 5) 



(A.16) 



= exp 



StiTi (x) exp 



P / d 3 x 



(VTTi 



+ 



J d 3 X^{x)J^{2 



J d 3 xd 3 x'J%(x)B~ 1 (x,x')J«(x') 



(A.17) 



/ d 3 x[J:(x)] 2 
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where B{x,x') = f38{x — x'). Using the expression (A. 2), (A. 3), (4.11), (4.12) the exponent 
in (A. 17) becomes: 

J d 3 x[J^x)] 2 

= J d 3 x i J dtid 3 xi J<f,(xi, ti)a n (x\ — x, t) i J dt2d 3 X2J ( f,(x2,t2)a 7r (x2 — x,t) 



Jdx ld X2Mx,)M^ l f££- [Ep t ~ ^ Sm [Ep tl ~ U)] e^(*-*) 



(27T) 



El 

v 



El 

v 



Inserting (A.18) and (A.16) in (A.4) and (A. 5) we obtain: 



(A.18) 



MM = exp \-\ j d^x'J,{x) { ^il^^ 1 ^^ 370 } J *V)\ ■ (A-19) 



Via the identity 



J dp 0_L e -i P °(t-t>) S{p 2_ m 2 ) 



cos [Ep(t - t') 

El 

v 



(A.20) 



we can rewrite Z^[J^\ as 



Zp[J<l>] = exp 



X - J d^x'J^A^x-x')^') 



(A.21) 



where 



A, 



^ 27T Mr? - m 2 \ P -iP-z 

(2vr) 4 p\p° 



This is what we have used in (4.15). 



B. Fields-Correlations Via Super-Field Projections. 



(A.22) 



In this appendix we will present some of the detailed calculations that we skipped in the 
text of the paper. Here, in particular, we will report the calculations which support the 
results in (5.29) and (5.30). Let us start from (5.29): 



J d8 1 d8 1 d8 2 d8 2 (0i0i) (8282) 



(2!) 2 zi 





Z2 



Using the fact that Q(z, z) = and the rules of Grassmannian integrations, we get for the 
expression above: 
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/ 



de 1 de 1 de 2 de 2 



(Mi) (M 2 ) J < 



d^xdOdO / d A x'd9'd9' 



A p (x! -x) + G [ n\xi -x)99 



+G { x\x - xi)Mi + G i p(x 1 - x)e l 6 + G { p{x - xi)^ij [A fi (x - x)] [A p (x - x') 
+Gp\x - x')9'9' + V - x)99 + G { p{x - x')99' + Gp\x' - x)9'9~\ A p {x' - x') 
A^x' - x 2 ) + Gp\x' - x 2 )9 2 9 2 + G ( p{x 2 - x')9'9' + G% \x' - x 2 )9'9 2 + Gp\x 2 - x')9 2 9' 
= \ j d9 1 d9 1 d9 2 d9 2 (Mi) (9 2 9 2 ) J d A xd9d9 J d A x'd9'd9' [Ap( Xl - x)Ap(x' - x 2 ) 

+A /3 (x i - x)G { p{x 2 - x')9'9' + G { p{ Xl - x)90Ap(x' - x 2 ) + G { *\ Xl - x)99 

Gp\x 2 - x')0'0'] Ap(x -x)x \Af,(x - x') + G ( p(x - x')9'9' + G ( p(x' - x)99 

+Gp\x - x')99' + Gp\x' - x)9'9\ Ap(x' - x') 

= \ j d9 1 d9 l d9 2 d9 2 (Mi) (M2) J d A xd9d9 J d A x'd9'd9' 

{A p { Xl - x)A p {x' - x 2 )Ap{x - x) [Ap(x - x') + G(z, z')\ Ap(x' - x') 

+ \A fi ( Xl - x)G ( p(x 2 - x')9'9'G { p(x' - x)99 + G { p( Xl - x)99Ap(x' - x 2 )G ( p(x - x')9'9' 

+G { p( Xl - x)99G ( p(x 2 - x')9'9'A p {x - x') A p {x - x)A p {x' - x')} 

= \ J d A x J d 4 x'Ap( Xl - x)A p {x - x)G ( p(x' - x)A l3 (x' - x')Gp ] {x 2 - x') 

+ \ J d A x J d A x'A l3 {x 2 -x)Ap{x-x)G [ p{x' -x)A p {x' -x^GPixY-x') 

+ \ J d A x J d A x'Gp\x 1 -x)Ap{x-x)Ap{x-x')A p {x' -x')gP{x 2 -x') 

(B.l) 

The three diagrams presented in the last expression are precisely 



(2!) 2 X! 



-• + 



x 2 (2!)2 Xl 



-• + 

X2 



(2!) 2 X! 



X2 



which is exactly the result we claimed in (5.29). 

Let us now give the proof of (5.30). Again we will make use of Q(z, z) = and of the 
standard rules of integration over Grassmannian coordinates. 
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J d9 1 d9 1 d9 2 d9 2 (Mi) (e 2 e 2 ) 



(2!) 2 zi 




(B.2) 



-J de 1 de 1 de 2 de 2 (Mi) (e 2 e 2 ) J d 4 xd9d9 J d 4 x'de'de' 



Ap( Xl -x) + G i ^ ) (x 1 - x)99 + G^\x - £i)Mi + G ( p( Xl - x)0 1 + G { £ \x - Xl )99~i 



x [A fi (x - x') + G(z, z')] [Ap(x' - x') + g(z', z')\ [Ap(x' - x) + Q{z\ z)] 



A^x - x 2 ) + G%\x - x 2 )9 2 9 2 + G^>(x 2 - x)99 + G K ^\x - x 2 )99 2 + G [ £>(x 2 - x)9 2 9 



= l -j d9 1 d9 1 d9 2 d9 2 (Mi) (9 2 9 2 ) j d 4 xd9d9 j d 4 x'd9'd9' 



ApQd - x)Ap{x - x 2 ) + Ap{ Xl - x)G^\x 2 - x)69 + G^\ Xl - x)99A p {x - x 2 ) 



+Gf } (xi - x)99G^>{x 2 - x)69\ [Ap{x - x') + g(z, z')] Ap(x' - x') [A p {x' - x) + G{z', z)} 

= l -j d9 1 d9 1 d9 2 d9 2 (Mi) (9 2 9 2 ) J d 4 xd9d9 J d 4 x'd9'd9' 

| [Ap( Xl - x)G^\x 2 - x)6d\ \a p {x - x')G^\x - x')9'9'j A p {x' - x') 

+ G ( ^\x 1 -x)99A p {x-x 2 ) 2Ap{x-x')G { k\x-x')9'9' Ap{x' - x')} 

= -J d 4 xd 4 x'A l3 (x 1 - x)G^\x - x')A p {x' - x')A p {x' - x)G^ ] {x 2 - x) 
d 4 xd 4 x'G^\x 1 - x)G^\x - x')Ap(x' - x')A p {x' - x)A p {x - x 2 ). 

(B.3) 



1 

2! x*T 



+ 



x 2 



2! 



Xl 



x 2 



This proves (5.30). 

In both of these proofs we have been rather pedantic with the details of the calculations. 
We did that because we wanted to help the reader not familiar with our formalism in going 
through the details of the calculations. 



C. Loop Identities. 

In this appendix we will prove various identities which, among other things, will help us in 
"proving" Eq. (5.48). We have put quotation marks around "proving", because ours will 
not be an analytic proof. Anyhow the identities that we will get will turn out to be rather 
usefull. Let us use the following notatio: dfi = d9d9, dfi± = d9\d9\, etc. The first identity 
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we want to prove is: 



J djxdm . . . dfj, n G(z, z 1 )C(z 1 ,z 2 ) . . . G(z n , z) = 0. 



(C.l) 



We shall call it the "first loop identity", because the strings of G makes a super- loop. Let 
us first observe that 



dzG(z 1 ,z)G(z,z 2 ) = I <fx G^ ) {x l -x)A^{x-x 2 ) + G^'(x 2 -x)A^(x 1 



/' 



d 4 x 


G%\ X1 


-x)G^\x 


- x 2 ) 


2 9 2 


d A x 


G^(x 2 


-x)Gf{x 


- Xl) 


Mi 


d A x 




-x)G%>{x 


- x 2 ) 


6*1^2 


d 4 x 


G { P(x 2 


-x)G^\x 


- Xl) 


6 2 6i 



If we add to the l.h.s. of Eq. (C.2) one more G and build the expression 



(C.2) 



J dzdz'G(z 1 ,z)G(z,z')G(z' ,z 2 ), 



(C.3) 



we can then do the integration over the Grassmannian variables and easily prove that: 



J dzdz'G(z 1 ,z)G(z,z')G(z',z 2 ) 

= J d 4 x [g^\ Xi - x)Ap(x - x')G^\x 2 - x') + A^(xi - x)G^\x' - x)G^\x' - x 2 ) 



+G^\ Xl - x)G^\x - x')A p {x' - x 2 ) 



I 



+ / d 4 x 



G^\x 1 -x)G^\x-x')G^\x'-x 2 ) 



2&2 



+ / d 4 x 



G < i\x 2 - x')G%>(x' - x)G^ ) (x - Xl 

± e 2 + / d 4 x 



xG { p(x'-x 2 ) 



71(71 



+ / d 4 x 



G^\x 2 -x')G^\x'-x)G^\x-x 1 ) 



G^>(x l -x)G^>(x-x') 
Q 2 Q\. 



(C.4) 



In general, with a string of all G, we can prove by induction that the following relation 
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holds (where z = zq, z' = z n+ \): 

j dzidz 2 ■ ■ ■ dz n G(z, zi)G(zi, 22) • • • G(z„, 2/) 
= (G^ *...* Ap i>i+1 * . . . * G^f ^ (x , x n +i) 



+ (G^f } * ... * G^fM (l ,Vl) fi n+l^+l 



(C.5) 



+ (G<f } *...*G ( / ) )(i n+ i 1 ioM 
+ (GSf } * ... * G<f } ) (x , X n+1 )0 O 0n+l 
+ (g^ * ... * G^) (x n+1 ,x o )0„ + i0 o . 
The various symbols appearing above are defined as follows: 

(G^ * ... * G^f^ (xo,x n+ i) = y d A X! . . . d A x n G ( x\x -x 1 )G ( x\x 1 -x 2 ) . . . G ( ^\x n -x n+1 ). 

(C.6) 

(we call this "convolution" .) While the first expression on the r.h.s. of Eq. (C.5) is defined 



as: 



(G^ * . . . * A^j+i * . . . * G^f ^ (x , x n+1 ) = 

= J d A xi . . . d A x n G^ ) (x - xi)G^f ) (xi - x 2 ) . . . G^f ) (xj_ 2 - Xj_i)G^ (xj_i - Xj) 
xA^x; - x i+ i)G)j 7 (x i+2 - X;+i)G}j ; (x i+3 - x i+2 ) ■■■G) i (x n - x n+1 ) 

(C.7) 

(we call this "oriented convolution"). Eq. (C.5) can be proved by induction as follows. Let 
us take P(n) = (Proposition of order n) to be given by Eq. (C.5). We already know, from 
the beginning of this appendix, that P(l) holds true [see Eq. (C.2)]. Next, we need to show 
that P(n) => P(n + 1). Indeed, take z' = z n+ \ in Eq. (C.5), multiply it by G(z n +i,z') and 
integrate over z n+ \. Using the definitions in Eq. (C.6) and Eq. (C.7) the results follows in 
a straightforward manner. 

Let bus now go back to (C.5), take z = z' and integrate over z, using the definitions 
(C.6) and (C.7), we easily get: 



dzdz\ . . . dz n G(z, z\)G(zi, z 2 ) ■ ■ ■ G(z n , z) = 0. (C8) 
A second loop equality which can be proved is 

dudm . . . dfi n G m (z, Zl )G mi (z u z2) ■ ■ ■ G m " (z n , z) = 0. (C.9) 



This is a generalization of (C.8) which is recovered when m = m\ = . . . = m n = 1. 
Using the property 4) of Eq. (5.17) we get that the l.h.s. of Eq. (C.9) is equal to: 
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/ 



dfidfii ■ ■ ■ d(i n G m (z, Zl )G mi (z!,z 2 ) . . . G m "(z n , z) 



I 



dfidfii . . . dfi n A™(x — xi) + mA™ l (x — xi)Q(z, z\\ 



A^(xi - x 2 ) + m 1 A^- 1 (x 1 - x 2 )g(z 1 ,z 2 ) 



(CIO) 



x ... x 



A™" 1 (x n -x) + m n AJ n 1 (x n -x)G(z n , Z ) 



= mm l ... m n Af-\x - Xi)A^~\xi - x 2 ) . . . A^~\x n - x) 
x J dfidfii . . . dfi n [Q(z,zi)G(zi,z 2 ) . . .<3(z n ,z)] = 0. 

In the steps above we have made use of (C.8), and of the relations: 

J dfidfii . . . dfi n g(z,zi)g(zi,z 2 ) . . . g{zi-i,Zi)A r ^ i {x i - x i+ i)g(z i+1 ,z i+2 ) . . . G(z n ,z) = 

(C.ll) 

and 

j dfidfii...dfi n Af(z,zi)A^(zi,z 2 )...A^(z n ,z)=0. (C.12) 

Relations (C.ll) and (C.12) can be easily proved considering the mismatch in the number 
of integration variables with respect to the variables present in the argument. 

We should not forget that we have to prove relation (5.48). The f( n \z, z') will be any 
combination of two-point functions with any number of loops inserted in order to reach the 
order n-th in perturbation theory. 

Let us start from n = 1 which is: 



1 

2! zT 



z 2 



FW(z,z'). 



It is easy to prove (5.48) for this F^, in fact: 
/ d9id6 2 d6id0 2 — 

J Zl 



z 2 



(C.13) 



Xl 



x 2 



= - j dzG^\x - xi)A p (x - x)G { *\x - x 2 ) = 



(C.14) 



This result comes from the fact that the integration contains dOdO, but the integrant does 
not contain any 9,9. 



Let us pass to the ¥^ 2 \z, z') functions. We could choose , for example, the following 

one: 



/ 




Xl 




X2 



(C.15) 



= - J dzdz'G { * ] (x - xi) [G(z, z')] 3 ) {x' - x 2 ) 
= - j d 4 xd 4 x'G^ ] (x - xi)GSfV - x 2 ) 



d/j,dfi'G s (z,z') 



= -d 4 xd 4 x'G^\x - x^G^ix' - x 2 ) J dfidfi'G 2 (z,z')G(z',z) 

In the last step we have used property 1) of Eq. (5.14). The piece in the last equality of 
Eq. (C.15) contained in square brackets is zero because of the loop identity in Eq. (C.8). 
So, everything turns out to be zero. Another example of ¥^(z, z') function is 



/ 



d0id0 2 dM02 





Xl 



X2 



dzdz'G { ^\x - xi)G(z, z')G{z\ z)G{z' , z)G ( ^ ) (x - x 2 ) 



(C.16) 



- -d 4 xd 4 x'G { x\x - Xi )G { k\x' - x 2 ) j d^'G(z,z')G(z',z')G(z',z) 



Once again the piece above contained in square brackets is zero because of the loop identity 
ofEq. (C.8). 

We checked many two-point functions of the third order in perturbation theory and 
we got zero because of the loop identity (C.8). We do not have an analytical proof to all 
orders, but we believe, from the many examples we worked out, that (5.48) is true. 

D. Jacobi Fields Correlations. 



In this appendix , to all orders in perturbation theory, we will prove thart: 

(^(xi)C*(x 2 )) j 8 = 



(D.l) 
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to all orders in perturbation theory. The reader could think that in C the and S' never 
interact and this could be the real reason behind (D.l). But, actually, the interact with 
4> and so, trough the (j), they could interact among themselves. 
We know from (5.39) that 

J d9 l d9 2 dhd9 2 9 2 9i{®{x l )${x 2 ))p = {d /> {x 1 )t^(x 2 ))p. (D.2) 

The l.h.s at the n-th order of perturbation theory will have the general form: 

J d9 1 d9 2 d9 1 d9 2 j dzdz / G(z 1 ,z)¥ ( - n \z,z')G(z',z 2 ) (D.3) 

where F^ is a function <D(g n ) that we do not write down explicitly. For the first and 
second order its explicit form is 

FW = Ig{z,z)S(z -z') (D.4) 



F< 2 ) = 



W(z,z') 



lG(z,z)G(z,z')G(z',z') 



dz"G(z, z")G(z", z")G(z", z)5(z - z') 



(D.5) 



The explicit form of (D.3) turns out to be: 



(D.3) 



ie l de 2 dd 1 de 2 9 2 e 1 J d^xdodo j d A x'de'de' 

A^xi -x) + G^\ Xl - x)99 + G^\x - xi)Mi 



+G^ ' (xi - x)6\6 + Gft'(x - x^Oh FW (z, z') 

x Ap(x' - x 2 ) + G^\x f - x 2 )9 2 9 2 + G ( k\x 2 - x')9'9' 

+G { p(x> - x 2 )9'9 2 + G { p{x 2 - x')9 2 9' 

dzdz' j d9 1 d9 2 d9 l d9 2 9 2 9i G i p(x 1 -x)9 1 9 ¥ {n) {z,z') G { p \x 2 - x')9 2 9' 
dzdz'99'G i p(x 1 - x)W {n \z,z')G { p{x 2 - x') 



(D.6) 

The last integral above is zero. The reason is the following: F^ n ^ comes from the 
integration of products of various Gr. As each Gr comes with an even number of 9 (where 
9 indicates either 9 or 9) attached , also F^ ra ^ will have an even number of 9 because the 
integrations take away always an even number of 9. As F^ n ^ depends only on z and z 1 and 
has an even number of 9, the only possible contributions must have the following strings 
of 9: 
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1) 666'6' which gives zero because 9' is attached also to the final G R in the second 
equality in (D.6). 

2) 99 which gives zero because of the initial G R (the term 9\9 in G R gives zero because 
of the projector #i#2 and the same for #2). 

3) 9' 9' same as in 2) but because of the final G R . 

4) 99' same as in 3). 

5) 99' same as in 2) 

E. Third Order Results For The Two-point Correlation. 

In this appendix we will provide the third order perturbative calculations for the (4>{xi)(f){x2)) p. 
correlation using the super-field technique. We will list below the super-diagrams that are 
needed together with the corresponding symmetry factors. These we calculate using the 
algebraic rules given in [29]. The set of super-diagrams needed are 
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Also, by considering the first and fourth diagrams in eq. (E.l), we note that the diagrams 
obtained by the interchange (z\ <R> z-i) gives a contribution. So all together we end up with 
ten superdiagrams at the third order ( In order to get the {(fxj>) components of these graphs 
we have to integrate them with the projector 61616262 like we did in (5.21) 
The calculations are long but straightforward and the result is 
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Of course using the set of super-diagrams (E.l) and different projectors we could get 
also the (cpX) and (cc) correlations up to the third order. 



F. Fluctuation-Dissipation Theorem Via Ward Identities. 

In this appendix we will present a non-perturbative derivation of the fluctuation-dissipation 
theorem and of other relations that we proved perturbatively in the paper. All these 
relations are basically Ward identities associated to the various symmetries present in the 
CPI, i.e., BRS, BRS, SUSY, etc. Most of the calculations presented in this appendix were 
done long ago by Martin Reuter whom we wish to warmly thank for his generosity. 

In this appendix we will not work in the proper KvN formalism with an Hilbert 
space, etc, etc, but in the Liouville one where a probability density p(ip) is used. 

Let us suppose we give the initial probability density p(— oo) at time t = — oo and that 
we want to evaluate the average of a quantity O. The path- integral expression for this 
average is: 

(C)p(-oo) = J ®<p®\®c®cQp(<p(-oo),c(-oo))exp (i J dtC^j (F.l) 

Let us ask for which symmetry transformation (whose infinitesimal transformation we 
indicate with 5) and for which p we have 

Wp(-oo) = 0. (F.2) 

Let us assume that the integration measure S>(p^\SicSic is invariant under the transfor- 
mation 5 and that the action S = J Cdt changes only by a surface term 

5S = ea{+oo) - ea(-oo), (F.3) 



where e is defined by (p e = cp + Se. 
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Let us change the variables of integration in (F.l) going over to those obtained via a 
symmertry transformation. We obtain: 

(C)p(-oo) = J S><p ■ ■ ■ (O + SO) [p (-00) + Sp (-00)] exp (iS + iea(+oo) - iea(-oo)) 
= J @ip...{Op + 05p + pSO + Op [iea{+oo) - iea(-oo)]} e*f 

(F.4) 

The first term inside the integral fives (O)p(_ 00 ) which is equal to the LHS because we did 
only a change in integration variables and so these two terms cancell. In this manner we 
get for the third term in the integral the following expression: 

(50)p(-oo) = -J@<p...O {5p(-oo) + p(-oo)ie [a(+oo) - a(-oo)}} e lS . (F.5) 
For 5 let us use the transformation generated by the super-symmetric charges [11] 

Q H = c a {i\ a - (38 a H) (F.6) 

Q H = c a {i\ a + (38 b H ) (F.7) 



where /3 is a real parameter which, when is set to zero, gives the BRS, BRS [11] charges . 
Let us now proceed to find the surface terms associated to this symmetry transformations: 



5 Qh C = -H3j t eN 



(F.8) 



(F.1 



where N = c a d a H and Q = c a u ab i\b- Let us now calculate the first term inside the brackets 
on the r.h.s of Eq. (F.5) for the transformation generated by Qh 



5p(-oo) = (S(p a + 5c a -^p\ (t = -00) = ec a (-oo)p ((p(-oo), c(-oo)) . 

The surface term inside the brackets of Eq. (F.5) gives 

iea(+oo) - iea(-oo) = /3e (c a d a H) (+00) - /3e {c a d a H) (-00). 
Inserting all this in (F.5) we get: 

d 



(F.10) 



(F.ll) 



I 



@ip...Oe 



c a (-00) 



f3c a (-00) d a H(v(-oo)) 



dip a (—oo) 

xp (<p(-oo), c(-oo)) - J @(pOpec a (+oo)d a H(+oo)p(-ooy 
= - j ®^OeQ H pe iS -p j 9<p... ec a d a HOpe i8 

For (3 = the second term vanish and the first is zero only if 

Qp = 



iS 



(F.12) 



(F.13) 
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where Q = Qh\/3=o = 0. The solution of (F.13) is 

P = p(<p)8(c). (F.14) 

Let us do the analog for the transformation generated by Qh ■ Calculating the various 
pieces we get 

(6 Qh 0)~ p = -J@<p...O {eQ H p) (-oc)e i5 ~. (F.15) 
This implies that (5q h O)~ p = if 

QhP = 0. (F.16) 

The solution of this equation was calculated in the second of Ref. [11] and is the canonical 
distribution 

p(<p,c)=e-P H 5(c). (F.17) 

So the averages taken with this p are exactly the thermal averages we have calculated 
ail-along in the core of the paper. 

Let us now choose as operator O the following one 

O = (p a (t)c b (0) ,t>0. (F.18) 

Doing the Qh transformation of O we get: 

5 Qh O = -eu ad c d (t)c\0) + <p a (t)eoj bd [i\ d (0) + f3d d H(0)] . (F.19) 

It was proved in [33] [34] [24] that 

(c d (t)c\0)) = (F.20) 

if t > 0, where the average is taken with the canonical average (F.17). So we get from 
(F.19) 

(5 Qh O) = = ecu bd (<p a (t) [i\ d (0) + (3d d H(0)}) (F.21) 

which implies: 

(<p a (t)i\ d (0)) = -P& a (t)d d H(0)). (F.22) 
From the equation of motion we have 

d d H(0) = oj de ^ e (0). (F.23) 

So (F.22) will become 

(<p a (t)i\ d (0)) = -/3^ e (^(t)v3 e (0)). (F.24) 

Now in general for time-independent Hamiltonian we have, for two observables A(t) and 
B(t), that 

(A(t 1 )B(t 2 ))=F(t 1 -t 2 ) (F.25) 

so 

±(A( tl )B(t 2 )) = -JL(A( tl )B(t 2 )) (F.26) 



- 55 - 



or 



(A{ tl )B{t 2 )) = -{A(h)B(t 2 )). 
Applying this last relation to (F.24) we get: 

{<p a (t)i\ d (0)) = -f3oj de {<p a (t)<p e (0)). 



(F.27) 



(F.28) 



and this is the fluctuation- dissipation theorem. 

With respect to the easier proof presented in [16] our proof has the virtue that it can 
be seen as a Ward-identity of Susy and that the canonical distribution is not put by hand, 
but it is derived [11] as a super-symmetric invariant state. 

Let us now prove few more identities. Let us choose p(<f,c) = 5 2n (c(— oo)). This 
distribution is BRS invariant [11] so: 



(5brsO) = J %> • • ■ (5brsO) S 2n c(-^)e iS = 0. 
where the BRS transformation [11] is 



(F.29) 



5ip a 


= ec a 


Sc a 


= 


Sc a 


= ie\ a 


5X a 


= 0. 



Let us take O = ip a (t)c b (0), then 

SO = 

So (SO) = implies 



= -ie [ic a (t)c b (0)-^ a (t)X b (0)]. 
{ip a (t)X b (0))=i{c a (t)c b (0))- 



(F.30) 
(F.31) 



This relation tell us that the dash-full propagator is equal to the Fermion one and to all 
order in perturbation theory. This confirms the statement that we made ail-along in our 
paper that this identity is related to the Susy invariance (better say to the BRS invariance 
which is the susy with j3 = 0) . 

Let us now choose as operator 



then 



and from (SO) = we obtain: 



6 = x a (t)c b (o), 

SO = ieX a (t)X a (0), 
(X a (t)X a (0)) =0. 



(F.32) 

(F.33) 
(F.34) 



This is "only" similar to the one we proved in the paper perturbatively to all orders, because 
there we proved it with the weight p = e~@ H S n (c), while here we proved it only for f3 = 0, 
but we think the proof given here should go through also for (5 ^ 0. 
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Let us now choose as operator 

6 = c a (t)^ b (0), (F.35) 

then 

SO = -ec a (t)c b (0), (F.36) 

so we get 

(c a (t)c b (0)) = 0. (F.37) 

This relation for a = b = 1 it was proved to all orders of perturbation theory in the paper 
and for the canonical distribution. Here we have given a non-perturbative proof. 

We like to conclude this appendix by saying that both the simplification in the dia- 
grammatics that we got via the introduction of super-diagrams and the various identites 
we found that will further simplify the diagrammatics, have their roots in the various 
symmetries of the CPI. 

G. Perturbative Check of The Fluctuation-Dissipation Theorem. 

In this appendix we will prove the relations (6.11) and (6.12). 
Let us remember that 

~X2)=J ^ ~ ^~ lV - X (G- 1 ) 

where p.x = p°t — p.x and p 2 = (p ) 2 — p 2 . 

Let us take the derivative with respect to ti, like on the l.h.s of Eq. (6.11) 

Let us remember that: 

^ = *(pP) _ *(_pP) = e(pP), (g.3) 

where we have defined the new symbol Q(p°). Let us also notice that: 

5 ((p ) 2 - f - m 2 ) = ^-5 (p° - Ep) -5(p° + Ep) , (G.4) 

and so we get: 



2Ep 



e(p°)<5 ((p ) 2 - f - m 2 ) = JL« (p° - Ep) -5(p° + Ep) . (G.5) 
We can now rewrite (G.2) as 

d?p p° 



d-A«(xi - Xo) ~ -i- ! d P P e -iEp{t2-ti) ip{Si-S 2 ) 

Hl ^{ Xl x 2 ) - (27r)32 ^ e e 

d 3 p p° 
j2^2E~p 



-i- j rf V P° e iE j frl-t2) e -ip.(xi-X2) 



(G.6) 
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Let us now turn to the r.h.s. of Eq. (6.11) and let us remember that Gr solves the equation: 

(□ + m 2 ) Gr = -S(x) (G.7) 

with Gr(x) = for t < 0. Let us then choose the contour Cr of integration along the real 
axis in the complex p plane in such a manner that the two poles are moved below the real 

Iim p° 



Ep 



Re p c 



axis. 

For t > we close the contour in the lower half plane: 

"°° d A p e-^ x f d A p e~^ x 



Gr(x) 



l-oo (2vr) 4 (pD + %e f - E} Jc R {^Y (p° + ief - E} ' 
So by the residue theorem we get: 



Gr{x) = i 



d 3 p 1 



-ipx 



d?p 1 e" 



-ipx 



(2vr)3 2E p (pO + ie f _ £| " J _ oo ( 2 vr)3 2% (p0 + , e ) 2 _ £| 
and comparing with Eq. (G.6) we obtain: 

d 

G R {x) = G R (xi - x 2 ) = (3—A(s(x 1 - x 2 ). 
For t < we get that Gr(x) is zero, so: 

In this relation let us now turn x — > — x and we get: 



So for t < we have: 



|A,(x) = -^(-x) 



(G.8) 
(G.9) 

(G.10) 

(G.ll) 

(G.12) 
(G.13) 
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Combining (G.ll) with (G.13) we can write: 

^A p (x) = ^G R (x) - ^G R (-x) (G.14) 

which is exactly the (6.11) or the FDT theorem at the zero order that is what we wanted 
to prove. 

We will now move on to prove Eq. (6.12). Let us first note that from (G.14) we obtain: 
G R (x! -x) = /39( Xl - x) J-A^xi - x). (G.15) 

Ol\ 

The first term on the l.h.s. of (6.12) has the following analytical expression: 



\lk * ^ x=\wj d4x ^ {Xl ~ X)A ^ X ~ X)G * ){X2 ~ X) 

= \j [ G k )( - Xi ~ x) ~ G R )( - X ~ A ^ x ~ X ) G R ] ^ ~ x ) 

(G.16) 

where we have made use of the relation (G.14). 

The second term on the l.h.s. of (6.12) has the following analytical expression: 



1 d 

•- 



= \f d A x^-G^\ Xl - x)A p {x - x)Ap(x - x 2 ). (G.17) 



2 dt\ x\ X2 
Let us now take the time derivative of (G.15): 

Qf i G < k\x 1 -x) = f38(h - t) w ^f){x x -x) + pe{ti ~ *)^2 A M*i ~ x )- ( G - 18 ) 



Moreover it is easy to prove that 



5(t 1 -t)-£-A fl (x 1 -x) = 0, (G.19) 

Oti 



and therefore, (G.18) becomes 



Wi G ( n\xi -x) = - t)-^Mxi ~ x). (G.20) 



As a consequence the (G.17) is turned into: 



Id \ J 1 „ r ,a „, ,d 2 



- - - \fi [ d 4 x0(t 1 -t)^A /3 (x 1 -x)A /3 (x-x)A^(x-x 2 ). (G.21) 

2 Oil X\ X2 2 J Ol\ 

Performing now an integration by parts and disregarding surface terms Eq. (G.21) can be 
rewritten as 
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1 d 



2 dti x x 



X2 



j d4x ^ tl ~ *^ *l ( Xl ~ X ) A P( X ~ x )A j3 (x2 ~ X) 



— X ) 

(G.22) 



where Ap ti {xi — x) = dA^Xj — x)/dti. Using (G.19) we get: 



1 d 



2 dt\ x\ X2 
If on the r.h.s. of the equation above we use (G.14), we obtain: 



j <fx6(ti - t)Ap tl (x! - x)Ap(x - x)Ap t2 (x 2 -x). 

(G.23) 



1 d 



2 dti x± X2 
Combining (G.16) with (G.24) we get: 



\\ j d A xG^\ Xl -x)A fj {x-x) [g^\x 2 - x) - G^\x - x 2 ) 



(G.24) 



1 d 



2 dti Xl 



1 d 



X2 2 dt\ Xi X2 

= \~pj dixG n\^ ~ X )M X " x )G ( r\x - x 2 ) 
~\\ I d i xG^\x-xi)A^x-x)G^\x2-x) 



2f3xi 



x 2 2/3xi 



X2 

(G.25) 

This is exactly the relation (6.12) which is the FDT of the first order in perturbation theory. 

Long but tedious calculations allows us to derive the FDT also at the second order 
in perturbation theory which is diagrammatically represented by the the three formulas 
(6.13), (6.14) and (6.15). 

For the first order contribution to the FDT an alternative and simpler derivation than 
the one presented here can be given either using the momentum representation or the 
super-field formalism. We will present them below: 

Momentum space derivation: Eq. (6.11), which is the FDT at zero order, is: 



d i 

—A/3(xi -x 2 ) = - [-iG R (xi - x 2 ) + iG R (x 2 - x t )} . 



(G.26) 
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Going into momentum space we get: 

74„ 



a 



- 1 ^ e ~ lp - {xl ' x2) ^p) 



d A p 



J (2vr 



-ip.(xi—x 2 ) 



-iG R (j>)} 



(G.27) 



or equivalently: 
d 4 p 



I 



(2tt)< 



_ipO e -ip-(xi-x 2 ) 



d 4 p 



-ip.(x 1 ~x 2 ) r_ 



iG R (p)] + 



d 4 p 



-ip.(x 1 -x 2 ) 



iG A (p)} 



(G.28) 



P\J (2vr) 4 1 nun J (2vr) 4 

where Ga is the advanced propagator. Formally we shall write the relation above as: 



—ip°Ap(p) = — [-iGr(p) + iGa(p)} , 



which diagrammatically is: 



p 



P 



Xi 



X 2 P 



Xl 



X 2 P Xl 



X2 



The FDT at first order (6.12) is given by: 



(G.29) 



(G.30) 



2! dh 



-• + •- 

X 2 Xl 



X2 



(G.31) 



11 
2! p 



Xl X 2 

The l.h.s. of (G.31) in momentum space is: 

q 



Xl 



X2 



2l ("^ 



p 



P 



+ 



P 



P 



(G.32) 



Let us now make use of (G.29) in the following manner: let us replace the first leg of the 
first diagram in (G.32) using the relation (G.29) and the same for the last leg of the second 
diagram in (G.32). Therefore, we get: 
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(G.32) = - 





(G.33) 



and this is exactly the momentum space representation of (G.31). 

One can notice how much simpler is this derivation than the one in configuration space 
that we presented earlier. 

Super-Diagram Derivation: 

Let us remember that: 



<(£(xi)<Kx 2 )) 



(i) 



dv\dv 2 



Z 2 



(G.34) 



where we have used the notation dv; L = dOidOiOiOi (i 
l.h.s. stays for first order in perturbation theory. 
It is then easy to prove that: 



1, 2) and the upper index (1) on the 



Zl 



Z2 



Xi 



I 



dv\ 



— • + • 

X 2 X\ 

d 9 (0 0*= 



X2 



=• ) + (•= 

z 2 o6 2 d6 2 zi 



) 



z 2 
(G.35) 

The first equality above derives from the fact that in x\ only the field <f> survives the 
integration over dv\. At the point z 2 in principle there could be any of the fields (f>, A, c, 
c, but only the correlations {(pcf)) and (4>X) survive because the (</>c), (cpc) are zero as we 
proved before in the paper. The second equality in (G.35) is easy to prove: the first term 
is exactly equal to the first term of the first equality. In fact, in z 2 only <\> survives which 
is the only term which does not carry extra 9 2 and in Z\ only (j> survives because of the 
integration dv\. An analogous reasoning applies to the second term. 

Let us now go back to (G.34) and do the time derivative of both, the l.h.s. and r.h.s. 



_d_ 



{^(x 1 )(f>(x 2 )) 



(i) 



dv\dv 2 




(G.36) 
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where with the arrow 
the r.h.s. of (G.36) as 



we indicate the time derivative. We can formally write 



- / dvidv 2 •= 
2 J zi 



=• = - I duidv 2 
Z2 2 J \Zi 



z z 



Z2 



(G.37) 



Using the relation (G.35) on the first term of the r.h.s. above, we get: 
d 



— {c/ ) (x 1 )c/ ) (x2)} W =\J dv x dv 2 



(G.38) 



Z2 



Let us analyze the two terms in the square brackets on the r.h.s. of (G.38). 
The first term is equal to 



I! 



dv\dv 2 



d_d_ 

_d9de v ""zi 



) 



Z2 



= \ f dv X dV2^B{*= 

2 J zi 



d_d_l 
86 88 1 ; 



Z2 



(G.39) 



1 

2xi 



X2 



where we have used integration by parts over Grassmann variables to obtain the first 
equality and we have performed the integration to obtain the second one. 

Similarly, the second term in the square brackets on the r.h.s. of (G.38) gives 



dv\dv2 



= 2 j du 1 dv 2 ( 



2 x\ 




(G.40) 
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Replacing now in (G.39) and (G.40) the continuous line with the arrow (which stays for 
the time derivative) with the difference of the dash-full line, like it is indicated in the FDT 
at zero order (6.11) , we get that: 



(G.39) 



2/? X! 



X2 



2/3 Xl 



X2 



(G.41) 



(G.40) 



2/? xT 




- • - 

X2 



2p 



Xl 



X2 



(G.42) 



Summing up the two relations above we get (6.12) which is the FDT to first order in 
perturbation theory. 

Both, the momentum and the super-field techniques can be used to go to higher order 
in the FDT. The calculations turn out to be simpler than in standard space, but still rather 
tedious, so we will avoid to report them here. 

This last part of this appendix has been somehow a demonstration of the many things 
one can prove by " manipulating" the Feynman diagrams without resorting to their analyt- 
ical expressions. These manupulations are "analogous" to what was done in the original 
" diagramma^ paper [1] for Quantum Field Theory, so we hope the use we made of the 
same word " diagramma^ is forgiven. 
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